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ABSTRACT 

We present and discuss the results of the Herschel Gould Belt survey (RGBS) observations in a ~11 deg^ area of the Aquila molecular 
cloud complex at J ~ 260 pc, imaged with the SPIRE and PACS photometric cameras in parallel mode from 70 jim to 500 jim. Using 
the multi-scale, multi-wavelength source extraction method getsources, we identify a complete sample of starless dense cores and 
embedded (Class O-I) protostars in this region, and analyze their global properties and spatial distributions. We find a total of 651 
starless cores, ~ 60% + 10% of which are gravitationally bound prestellar cores, and they will likely form stars in the future. We also 
detect 58 protostellar cores. The core mass function (CMF) derived for the large population of prestellar cores is very similar in shape 
to the stellar initial mass function (IMF), confirming earlier findings on a much stronger statistical basis and supporting the view that 
there is a close physical link between the stellar IMF and the prestellar CMF. The global shift in mass scale observed between the 
CMF and the IMF is consistent with a typical star formation efficiency of ~40% at the level of an individual core. By comparing 
the numbers of starless cores in various density bins to the number of young stellar objects (YSOs), we estimate that the lifetime 
of prestellar cores is ~ 1 Myr, which is typically ~ 4 times longer than the core free-fall time, and that it decreases with average 
core density. We find a strong correlation between the spatial distribution of prestellar cores and the densest filaments observed in the 
Aquila complex. About 90% of the Herschel-identihcd prestellar cores are located above a background column density corresponding 
to Ay ~ 7, and ~ 75% of them lie within filamentary structures with supercritical masses per unit length > 16 Mo/pc. These findings 
support a picture wherein the cores making up the peak of the CMF (and probably responsible for the base of the IMF) result primarily 
from the gravitational fragmentation of marginally supercritical filaments. Given that filaments appear to dominate the mass budget 
of dense gas at Ay > 7, our findings also suggest that the physics of prestellar core formation within filaments is responsible for a 
characteristic “efficiency” SFR/M^ense ~ x 10”^ yr"^ for the star formation process in dense gas. 
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1, Introduction: The Herschel Gould Belt survey 


Understanding how dense cloud cores and protostars form out 
of the diffuse interstellar medium (ISM) is a fundamental ques¬ 
tion in contemporary astrophysics (e.g., McKee & Qstriker|2007| 
and other recent reviews in |Beuther et al.|2014 ). Much progress 
is being made on this front thanks to imaging surveys with the 
Herschel Space Observatory ( [Pilbratt et al.|2010|). Its f ar-infrared 
and sub millimeter cameras PACS ( [Poglitsch et al.||201()| and 
SPIRE ( [Griffin et al. 112010! ) provide unprecedented sensitivity 
and dynamic range at wavelengths around the peak of the spec¬ 
tral energy distributions (SEDs) of starless cores and protostars. 

In particular, the bulk of nearby (J 500 pc) molecular 
clouds, mostly located in the Gould Belt (e.g., |Guillout|[2001[ 
Perrot & Grenier||2003| ), have been imaged at five wavelengths 
between 70 pm and 500 pm as par t of the Herschel Gould Belt 
survey (RGBS - Andre et al. [20 1 0} . Observationally, the molec¬ 
ular clouds of the Gould Belt are the best laboratories at our 


^ Herschel is an ESA space observatory with science instruments pro¬ 
vided by European-led Principal Investigator consortia and with impor¬ 
tant participation from NASA. 


disposal for investigating the star formation process in detail, 
at least as far as low-mass stars are concerned. They are the 
only clouds for which the ~ 15" angular resolution of Her¬ 
schel around T ~ 200 pm is sufficient to resolve th e typical Jeans 
length ~ 0.03 pc in cluster-forming clumps (e.g.,|Larson|1985[ 
|My ers TT^ . 

The 15 or so nearby clouds covered by the RGBS span a 
wide range of physical and environmental conditions, from very 
active, cluster-forming complexes such as the Orion A & B gi¬ 
ant molecular clouds (GMCs) or the Aquila Rift cloud complex 
(e.g., [Dame et al.|200H|Gutermuth et al.|2008| ) to quiescent re¬ 
gions with no star formation activ ity whatsoever, such as the Po- 
laris flare translucent cloud (e.g., |Reithausen et aH2002 


Ward- 


[Thompson et al.|20T0| ). The total surface area covered by the sur- 
vey exceeds 160 deg^ (cf. http://gouldbelt-herschel.cea.fr/ for the 
list of all target regions). The RGBS will eventually provide an 
essentially complete census of (solar-type) prestellar cores and 
Class 0 protostars with well-characterized luminosity and mass 
functions in most nearby star-forming regions. 

The main scientific goals of the RGBS are to clarify the na¬ 
ture of the relationship between the prestellar core mass func- 
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tion (CMF) and the stellar initial mass function (IMF) and to 
elucidate the physical mechanisms responsible for the growth of 
structure in the cold interstellar medium (ISM), leading to the 
formation of prestellar cores and protostars in molecular clouds. 

Initial results from the RGBS have already been presented 
in several “first-look” papers and may be summarized as fol¬ 
lows. The RGBS observations confirm the omnipresence of fil¬ 
aments in nearby molecular clouds and suggest an intimate con¬ 
nection between the filamentary structure of the cold ISM and 
the formation process of prestellar cores ( [Andre et al.||201Q] 
Men’shchikov et al.|2010| ). While molecular clouds were already 


known to exhi bit large-scale filamentary structures long before 
Schneider & Elmegreen||197^ |Rartmann||2002t 


Herschel (e.g., 

Myers||2009j), the Herschel observations from the RGBS (e.g., 


Men’shchikov et aL]|2010| [Miville-Deschenes et al.||2010|) and 


other imaging surveys such as RiGAL (jMolinari et al.pOlQ 
Schisano et al. ] |M4li, HOBYS dMotte et al.||201Ul IHilleTal 


2011 ), and EPoS ( Henning et al. 2010| l now demonstrate that 


these filaments are truly ubiquitous in the cold ISM, present a 
high degree of universality (e.g., [Arzoumanian et al.pOlT]), and 


likely play a central role in the star formation process (see [Andre 
. [2014[ for a recent review). In any given cloud, Herschel 


[et al 


[201 l[[Palmeirim et al.[2013t[Alves de Oliveira et al.[2014[ ). Other 

groups have reported results in broad agreement with ou r RGBS 
finding of a common filament width. [Juvela et al.[ ( [2Q12[ ) found a 
typical FWRM width of ~ 0.2-0.3 pc for a number of filaments 
mapped as part of the Herschel “Galactic Cold Cores” project in 
clouds with (rather uncertain) distances ranging from ~ 100 pc 
to a few kpc. [Ysard et al.[ ( [2^013[ ) reported a mean FWRM width 
~0.1 pc for the LI506 filament in Taurus but found significant 
variations - by up to a factor of ~ 2 on either side of the mean 
width - along the length of the filament. [Smith et ah] ( [2014[ ) ex¬ 
plored filament properties in a set of numerical hydrodynamic 
simulations and found a range of filament widths rather than a 
constan t value. Recent magneto -hydrodynamic (MRD) simula¬ 
tions by [Ntormousi et al.| ( [2015[ ), however, suggest that non-ideal 
MRD turbulence can account for the properties of observed fila- 
ments much bette r than hydrodynamic turbulence does (see also 
[Rennebelle|2013[ ). 

The origin of the common inner width of interstellar fil¬ 
aments in nearby clouds is not yet well understood. A pos¬ 
sible interpretation is that filaments result from planar inter¬ 
secting shock waves due to supersonic interstellar turbulence 
(e.g. [Pudritz & Kevlahan|[2013[ ). and that the filament width 
corresponds to the sonic scale below which the turbulence be¬ 
comes subsonic in diffuse, non-star-forming molecular gas (cf. 
[Padoan et al.[[200l| ). Alternatively, a characteristic width may 
arise if interstellar filaments are formed as quasi-equilibrium 
structures in pressure balanc e with a typical ambient IS M pres¬ 


sure Pext^^-SxlO"^ Kcm ^ (Fischera & Martin 2012 S. Inut- 


suka, private communication). Yet another possibility is that the 
filament inner width may be set by the dissipation mechanism of 
MRD waves due to ion-neutral friction ( [Rennebelle|2Q13[ ). 

The early results from the RGBS further suggest that prestel¬ 
lar c ores and protostars form primarily in the de nsest filaments 
(e.g., [Andre et al.[[2010t [Polychroni et aL][2013[ ), for which the 
mass per unit length exceeds the critical line mass of nearly 
isothermal, long cylinders (e.g., [Inutsuka & Miyama|[1997[ ), 
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Mine,crit = 2cJ/G ~ 16Mo/pc, where Cs ~ 0.2 km/s is the 
isothermal sound speed for molecular gas at T ~ 10 K. They 
also confirm the existence of a close relationship between the 
prestellar CMF and the stellar IMF in the regime of low to in¬ 
termediate stellar masses (~0.1-5Mo - [Konyves et al.[[2010[ ). 
These Herschel findings support a scenario according to which 
the formati on of solar-type stars occurs in two main steps ( [Andre[ 
et al.|2014| ): first, the dissipation of kinetic energy in large-scale 


magneto-hydrodynamic (MRD) flows (turbulent or not) gener¬ 
ates a quasi-universal web-like filamentary structure in the ISM; 
second, the densest filaments fragment into prestellar cores (and 
ultimately protostars) by gravitational instability. 

In this paper, we present the “first-generation” catalog of 
dense cores obtained from RGBS data in the Aquila Rift cloud 
complex and discuss the global properties of these dense cores 
in relation to the filamentary structure of the complex. In par¬ 
ticular, we use these results to quantify the role of filaments in 
the star formation process. The present study extends and re- 
inforces our early Herschel findings in Aquila ( [Konyves et al 


imaging reveals a whole network of filaments, and a detailed 
analysis of the radial column density profiles of the nearby, re¬ 
solved filaments observed in the RGBS shows that they are char¬ 
acterized by a very narrow distribution of central widths with a 
typical full width at half maximum (FWRM) value ~ 0.1 pc 
and a dispersion of less than a factor of 2 ([Arzoumanian et al] 


[201Q[ [Andre et~ar][201Q[ ) on the basis of a more advanced ex¬ 
amination of the data with improved data reduction, source ex¬ 
traction, and source characterization. The paper is organized as 
follows. Section introduces the Aquila Rift region. Section 
provides details about the Herschel imaging observations and 
the data reduction. Section [^presents the dust temperature and 
column density maps derived from Herschel data, describes the 
filamentary structure seen in these maps, and explains how dense 
cores were extracted, prestellar cores selected, and their proper¬ 
ties measured from the maps. In Sect. we discuss estimates 
of prestellar core lifetimes, the observational evidence of a col¬ 
umn density threshold for prestellar core formation, the spatial 
distribution of extracted dense cores, and the strong connection 
with the filamentary structure of the Aquila cloud. We also com¬ 
pare the CMF of the Aquila sample of prestellar cores with the 
IMF, and link the global star formation rate of the complex with 
the total mass of dense gas above the column density threshold 
for star formation. Finally, Sect, [^concludes the paper by sum¬ 
marizing the RGBS results in the Aquila region and discussing 
possible implications for our understanding of star formation on 
GMC scales. 


2. The Aquila Rift region 

The Aquila Rift molecular cloud complex corresponds to a large 
extinction feature (see Prato et al. 2008), located above the 
Galactic plane (Z? ^ 4°) at galactic longitudes between I = 30° 
and I = 50°. The portion of the cloud complex mapped with 
Herschel as part of the RGBS corresponds to t he western high- 
extinction area of the Aquila Rift at I < 35° (see[Bontemps et al. 
[MTO] ). 

While the northern part of the Aquila high-extinction area 
harbors the well-documented Serpens Main star-forming region, 
the properties of the southern part (the focus of the present paper) 
remained largely unexplored until Spitzer infrared observations 
(e.g.,[Gutermuth et al.[[2008| ). This extinction-defined area (see 
[Bontemps et al.|2010[ ), rich in gas but initially thought to be al¬ 
most devoid of star formation ( [Prato et al.|2008[), is now k nown 
to harbor two cluster-forming clumps ( Maury et al.|2011[ ): Ser¬ 
pens South, a young protostellar cluster showing very active re¬ 
cent star formation and embedded in a dense filamentary cloud 


([Gutermuth et al. [[2008[[Bontemps et aL][2010[[Nakamura et al. 


2011 


2013 


jTeixeira et al.[[26l2t [Friesen et ai.[[2013| [Kirk, R. et al. 

[Tanaka et al.|2013[ )7^and W40 a young star cluster associ¬ 


ated with the eponymous RII region, also known as Sharpless 2- 
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64 dSmith et al.||1985l |Vallee||1987l |Kuhn et al.||2010l |Pirogov| 
let al.|2013H . 

Whether or not the southern part of the Aquila high- 
extinction region and the Serpens Main cloud are at the same PACS data reduction 
distance is still a matter of debate ([Bontemps et al.|20TQt [Maury 


ping took ~12 hours. The above strategy is similar for all the 
parallel-mode SPIRE/PACS observations of the RGBS. 


et al. 


et al. 


2011[ Loinard|2013| ). Based on stellar photonietry, |Straizy's 


(2003) concluded that the front edge of the Aquila molecu¬ 


lar cloud was at 255+55 pc. Using VLBI observations, however, 
Dzib et al^ ( |2010| ) measured the trigonometric parallax of the bi¬ 


nary system EC95 in the Serpens Main region and obtained a 
distance of 415+15 pc. Erom the extinction maps obtained by 
[Bontemps et al.| ( [20T0| , the respective extinction features toward 
the eastern Aquila Rift region (containing Serpens South, W40, 
and Sh2-62) and the Serpens Main cloud are seen as clearly 
distinct regions. It is therefore possible that the two clouds are 
not physically associated, but located along neighboring lines 
of sight. While the method used by |Straizys et al.|p003||1996) 
would naturally be sensitive to the first dust extinction screen 
along the line of sight, the larger VLBI-based distance of the 
Serpens Main core by |Dzib et ah] ( |201Q| ) suggests that Serpens 
Main is located behind the extinction wall associated with the 
Aquila clouds (Serpens South being the highest extinction region 
found inside the Aquila Rift complex). A distance of 260 pc for 
the Aquila Rift complex also suits the MWC297/Sh2-62 region 
since the young star MWC297 itself has an accepted distance 
of 250 pc ( [Drew et ar]|1997| ). It is finally worth noting that the 
visual extinction map derived by |Cambresy| ( |1999| ) from opti¬ 
cal star counts and only tracing the first layer of the extinction 
wall has exactly the same global aspect as the 2MASS extinc¬ 
tion map of |Bontemps et al.| ( |2010| ), suggesting that both Serpens 
South and the W40/Aquila rift/MWC297 region are associated 
with this extinction wall at 260 pc. We will thus adopt a dis¬ 
tance d = 260 pc for the entire Aquila complex, throughout this 
paper (see Appendix however, for a brief discussion of how 
our results would change had we adopted a distance of 415 pc 
instead). 

Rescaled to our adopted distance, the molecular mass of the 
entire Aquila Rift has been estimated from CO an d extinction 


data to be 2-5X10^ Mq wit hin a 25 pc-radius region ( Dame et al 
1987[ [Straizys et al.||2003] ). Rescaled to the same distance, the 


virial mass for the entire Aquila Rift estimated by [Dame & Thad-| 
deus ( 198^ is ~ 3.3x10^ Mo, suggesting that the whole c omplex 


is close to virial balance on large scales. More recently, [Tanakal 
et al.|(|2013]) obtained a virial parameter ~0.08-0.24 for the Ser¬ 


pens South filament (ag ain rescaled to a dista nce of 260 pc) on 
~ 0.5 p c scales (see also [Kirk, H. et al.|2013| ), and [Maury et al.| 
( 201 1| derived a high star formation rate of ~23 M© Myr“^pc“^ 
for the protocluster associated with the filament (of total mass 
~610 Mq, also using d = 260 pc). Altogether, these results 
suggest that the Aquila Rift complex is globally gravitationally 
bound on scales of ~25 pc and includes a few highly unstable 
(sub-virial) clumps on the verge of forming rich star clusters on 
sub-parsec scales. 

3. Observations and data reduction 

The Herschel Gould Belt survey observations of the Aquila 
Rift complex were taken on 24 October 2009 during the Sci¬ 
ence Demonstration Phase of Herschel ( [Pilbratt et al.[[2010| . 
The SPIRE/PACS parallel-mode scan maps co vered a common 
~11 deg^ area both by SPIRE ( [Griffin et al.||2010| and PACS 
(|Poglitsch et al.||2010|(. With one repetition in two orthogonal 
observing directions (OBSIDs: 1342186277, 1342186278), the 
scanning speed was 60"s“*, and the total duration of the map¬ 


The individual scan directions of the parallel-mode PACS data at 
70 yum and 160 /im were reduced with RIPE ( [Qtt|201l| version 
9.0.3063, provided by the Herschel Science Center. 

Starting from the raw data (level-0) and up to the level-1 
stage, standard steps of the default pipeline were applied. The 
PACS photometer flux calibration scheme was applied using the 
up-to-date r^onsivity and correction factors (PACS ICC report, 
Balog et al.|Jof the executed RIPE version with the calibration 
file set PACS_CAL_45_0. During the actual processing of the 
data, we created masks to avoid bad and saturated pixels, cali¬ 
bration blocks and their unexpected transient effect on the sub¬ 
sequent frames. Besides the flat-field correction, we applied a 
non-linearity correction to the data (PACS ICC report. Billot et 
al.)^. The PACS bolometers enter a non-linear regime for point 
sources above ~100 Jy/beam in all bands (70/100/160 yum), and 
the flux densities of brighter targets are underestimated by typ¬ 
ically a few percent. The applied non-linearity correction of the 
PACS bolometer signal had a very minor effect on the Aquila 
data. Cosmic ray hits on the detectors were removed with the 
"second-level deglitching" method of RIPE. To make best use 
of the deglitcher, we took special care to prepare its input data. 
Eirst, a high-pass filtering with a scan-leg length outside of a 
protective object mask was performed. Next, the second-level 
deglitching was then applied on these temporary data. Baseline 
subtraction was only used for deglitching purposes, but not on 
the resulting level-1 frames. The slew/turn-around data at the end 
of the scan legs were also preserved in the processing. 

Eurther treatment of the flux- and pointing-calibrated level- 
1 time series and the projection of the combined scans were 
perform ed with an IDE -based map-maker, Scanamorphos, ver¬ 
sion 20 ( Roussel 2013| p The processing is fully automated with 
some user-defined options. It consists of the main functionalities 
of subtracting both the thermal and non-thermal components of 
the brightness drifts, as well as detecting and masking remain¬ 
ing glitches and brightness discontinuities in the PACS data. In 
the final map projection, we adopted a spatial grid of 3'7pixel. 
Scanamorphos also provides associated maps of error, total drift, 
and weight (see Sect. 3.7 of [Roussel|2013[ for details). The error 
map provides the error on the mean brightness in each pixel. 
In the case of two scan directions, an additional “clean” map is 
produced, which is a signal map weighted so that noisy scans are 
excluded for each pixel. The clean map is only used for diagnos¬ 
tic purposes. In the PACS map processing, a final step was per¬ 
formed to remove long artifact glitches, which remained mainly 
in the 70 yum map, due to a jump in the brightness of the PACS 
data that could affect whole array rows. Thanks to the various 
planes of the output map, we could replace only affected pixels 
by “clean map” pixels. Our PACS output (level-2) fits files were 
produced in Jy/3"-pixel units. 

Eor PACS data, the absolute flux accuracy of point sources is 
3% in the blue band (70 yum) and better than 5% in the red band 
(160 jam) (cf. PACS ICC report by Muller et al.^). The extended 
source calibration is more uncertain. In this paper, we conserva¬ 
tively adopted 10% and 20% absolute calibration uncertainties 
for the integrated source flux densities measured in the 70 jam 
and 160 yum bands, respectively (see also Sect.[4.6[below). 


^ http://herschel.esac.esa.int/twiki/bin/view/Public/PacsCalibrationWeb 
^ The documentation and repository of the software can be found at: 
http://www2.iap.fr/users/roussel/herschel 
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SPIRE data reduction 


The SPIRE 250 //m, 350 /iin, and 500 pm data were reduced 
with HIPE version 10.0.2751 using modified pipeline scripts. 
The nominal and orthogonal scan directions were processed in¬ 
dividually, and combined in a second step. Data taken during the 
turnarounds of the satellite were not included in the final maps. 

The raw level-0 data (in engineering units) were processed 
to level-0.5 (in physical units) using the relevant calibration trees 
(SPIRE_CAL_10_1) built in HIPE. The following pipeline steps 
to level-1 (cf. [Dowell et al.|20T0| consist of: 1) converting detec¬ 
tor timelines to angles on the sky, 2) creating the pointing prod¬ 
uct for the observation, 3) correcting for thermistor-bolometer 
electrical crosstalk, 4) correcting temperature drifts and detect¬ 
ing temperature jumps, 5) identifying glitches caused by cosmic 
rays, for which the assumption was that all glitches affect all 
bolometers of SPIRE simulataneously, 6) applying the low-pass 
filter response correction, 7) applying the flux conversion, and 
8) searching and correcting for cooler burps by recalculating the 
temperature drift calibration table. (A cooler burp is a steep tem¬ 
perature rise which reaches a stable plateau ~6-7 hours after the 
cooler recycle ends.) 

As the Aquila region is dominated by extended emission 
from the interstellar medium, relative gain factors appropriate 
to extended sources were applied to the bolometer timelines. 
These gains, determined by the SPIRE ICC, represent the ratio 
between the response of each bolometer to the extended emis¬ 
sion and the average response. Variations in the specific response 
of each bolometer arise due to variations in the beam area among 
bolometers. 

The destriper module of the pipeline was used in an iterative 
manner. The iterative process starts with level-1 timelines for 
both scan directions, and reconstructs an initial naive map which 
is only corrected for a median offset. The destriper then fits a 
constant level to the difference between each input timeline and 
the corresponding map timeline, subtracts the fit from the orig¬ 
inal timeline, and reconstructs another map. By default, bright 
sources are excluded during baseline fitting. These steps are iter¬ 
ated until convergence. We adopted default grid pixel sizes of 6", 
10", 14" for the SPIRE 250 yum, 350 yum, 500 pm wavelengths, 
respectively. The output (level-2) fits files for each SPIRE wave¬ 
length were in Jy/beam units. Eor SPIRE data, the absolute flux 
accuracy is better than ~ 5% for point sources ( jBendo et al 
2013|) and better than ~ 10% for extended sources (cf. [Griffin 


et al. 


2013| ) in the three bands. 


Map-making tests and consistency of the SPIRE PACS maps 

SPIRE and PACS map-making tests and benchmarks were car¬ 
ried out in early 2012 by SPIRE/PACS ICC members, map- 
maker developers, and Herschel Key Program representatives. 
The public SPIR^ and PAC5|^ results of this test campaign, 
which compared the performance of several publicly available 
map-making methods, justify our choice of the destriper pipeline 
with 0^^-order baseline removal (PO) for SPIRE data reduction 
and the choice of Scanamorphos for PACS map-making. In par¬ 
ticular, the destriper PO, the default map-maker in the SPIRE 
scan-map pipeline since HIPE v9, performed remarkably well 
and compared favorably among all map makers in all test cases 
except for those suffering from the “cooler burp” effect. Eur- 
thermore, the destriper can handle observations with complex 


^ https://nhscsci.ipac.caltech.edu/sc/index.php/Spire/SPIREMap- 
MakingTest2013 

http://herschel.esac.esa.int/twiki/bin/view/Public/PacsCalibrationWeb 
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extended emission structures and with large-scale background 
gradients very well. Power-spectrum tests carried out on SPIRE 
scan maps by the SPIRE ICC (see also [Miville-Deschenes et al. 
[2010[ for the case of the HGBS images of the Polaris flare cirrus 
cloud) demonstrate that large SPIRE maps such as the HGBS 
maps trace a wide range of angular scales reliably, from ^ 30' 
or more down to the SPIRE angular resolution (e.g. ~ 18" at 
250 pm). This high spatial dynamic range is a key advantage 
of the Herschel/SFIRE images (compared to, e.g., ground-based 
submillimeter continuum data), which makes our Herschel sur¬ 
vey simultaneously sensitive to both large-scale structures in 
molecular clouds (e.g. filaments) and small-scale structures such 
as individual prestellar and protostellar cores. 

As for PACS maps, comparison metrics showed that the pho¬ 
tometry of both point-like and extended sources carried out on 
Scanamorphos maps is highly consistent with the results ob¬ 
tained on maps produced with other map-makers. 

The relative astrometry between the SPIRE and PACS im¬ 
ages was tested by cross-correlating the various maps after re¬ 
projecting them on the same grid. Using the evolved pointing 
products in the Herschel system, a good match was found be¬ 
tween the intensity peaks seen in the SPIRE and PACS maps on 
a resampled 2"/pixel scale grid. 

The 3 "/pixel maps used in the present analysis are thus well 
registered and have a relative astrometric consistency better than 
2". The absolute astrometry of the Herschel images was also 
compared with publicly-available Spitzer data, as well as high- 
positional accuracy (<!") 3 mm IRAM Plateau de Bure observa¬ 
tions of a small field at the center of the Serpens South filament 
( [Maury et al.[[2011[ ). The final absolute astrometric accuracy of 
the Herschel maps is estimated to be better than 3". 

The parallel-mode PACS and SPIRE maps used in this paper 
were all converted to MJy/sr units and reprojected to a common 
3" pixel grid. The conversion of the PACS maps from Jy/3"- 
pixel units to MJy/sr units was obtained using a square pixel 
area of 9 arcsec^. Eor the SPIRE unit conversion from Jy/beam 
to MJy/sr units, we assumed the beam areas measured in 1"- 
pixel beam maps by the SPIRE ICC, as given in Table 5.2 of 
the SPIRE Observer’s Manual v.2.2 (29 Nov. 2010), namely 
426 arcsec^, 771 arcsec^, 1626 arcsec^ at 250 pm, 350 pm, 
500 pm, respectively. The half-power beam width (HPBW) res¬ 
olutions of the maps are 8.4", 13.5", 18.2", 24.9", and 36.3" 
at IQ pm, \6Qpm, 25Qpm, 350yum, and 500yum, respectively. 
These high-quality maps are publicly available from the Her¬ 
schel Gould Belt Survey Archiv^ 

4. Results and analysis 

4.1. Dust temperature and column density maps 

We used the Herschel images to construct an H 2 column density 
map (/Vh 2 5 Fig- ^ temperature map (T^, Eig. of 

the Aquila field. We first smoothed all Herschel images (repro¬ 
jected to the same 3" pixel grid - see above) to the 36.3" HPBW 
resolution of the SPIRE 500 pm data. 

A zero-level offset, obtained by correlating the Herschel data 
with Planck and IRAS data (cf. [Bernard et al.[[2010| ), was also 
added at this stage to each Herschel map. The added offset val¬ 
ues were 27.7, 159.8, 169.6, 94.4, and 41.4 MJy/sr at 70, 160, 
250, 350, and 500 pm, respectively. Assuming optically thin dust 
emission at a single temperature for each map pixel, we then 
fitted a modified blackbody function of the form f = By(Tf)KyL 

^ http://gouldbelt-herschel.cea.fr/archives 
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Fig. 1. a) H 2 column density map of the Aquila region at 18.2" angular resolution, as derived from RGBS data using the method described in 
Sect. |4.1| b) Same map as in the l eft pa nel with the p ositions of the 446 candidate prestellar cores and 58 protostellar cores identified in the Herschel 
images with getsources (see Sect. |4.4| and Sect. |4.5| ) shown as black and magenta triangles, respectively. Yellow triangles locate additional prestellar 
and protostellar cores which were excluded from the analysis and discussion of this paper, due to likely contamination by m ore distant objects 
belonging to background CO clouds at significantly higher LSR velocities than the bulk of the Aquila complex (cf. Sect. |4.3| ). The orientation of 
the galactic coordinate axes is indicated at the lower right of each panel. The lower left edge of the map is oriented almost parallel to the galactic 
longitude axis at Gb ~ 2° above the Galactic plane. 



Fig. 2. Dust temperature map of the Aquila region at 36.3"resolution, as derived from RGBS data (see Sect. |4.1| ). The white contour outlines the 
W40 RII region which has > 17 K. 


to the four observed data points from 160 jim to 500 fim on a 
pixel-by-pixel basis, where /y is the surface brightness at fre¬ 
quency y and By{T^) is the Planck blackbody function. Each SED 
data point was weighted by l/cr^ap where (Xcai corresponds to the 
absolute calibration error (20% of the intensity at 160 yum and 
10% for the SPIRE bands). We adopted a power law approxima¬ 
tion to the dust opacity law per unit mass (of dust-hgas) at sub¬ 


millimeter wavelengths, namely /c^ = 0.1 x (/l/300yum)“^ cm^g, 
and fixed the dust emissivity index yS to 2 (cf. |Hildebrand| 198^ ). 
Based on a detailed comparison of the Herschel results with the 
near-infrared extinction study of [Alves et al.| ( |200T] ) for the star¬ 
less core B68, [Roy et al.| ( |2014| ) concluded that these dust opacity 
assumptions are likely appropriate to better than 50% accuracy 
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over the whole range of H 2 column densities between ~ 3x10^^ 
cm“^ and ~ 10^^ cm“^. 


In the SED fitting procedure, the gas surface density distri¬ 
bution (E) and the dust temperature were left as two free param¬ 
eters. The H 2 column density (A^H 2 ) was then calculated from 
E = adopting a mean molecular weigh t per hydrogen 

molecule = 2.8 (e.g., Kaulfmann et al.|2008| ). Based on this 
SED-fitting method, we derived both a standard column density 
map at the ~ 36.3" resolution of the SPIRE 500 pm data and 
a ‘high-resolution’ column density map at the ~ 18.2" resolu¬ 
tion of the SPIRE 250 pm data. The procedure used to construct 
the ‘high-resolution’ column density map is based on a multi¬ 
scale decomposition of the imaging data and described in detail 
in Appendix A of |Palmeirim et al.| ( 2013| ). 


Both the standard and the high-resolution column density 
maps were tested against a near-infrared extinction map of the 
Aquila/Serpens region derived from 2MASS data (see |Bontemps| 
|et al.|2010 , [Schneider et al.|2011 ), the latter with a EWHM spa¬ 
tial resolution of -120". To do this, the Herschel column density 
maps were smoothed to 120" and converted to vis ual extinction 
units assuming (cm“^) = 0.94 x 10^^ Ay (mag) (Bohlin et al 


1978|). We then derived ratio maps of the converted Herschel 


maps to the Ay map from 2MASS on the same grid. In most 
of the field covered by Eig.[2 the ratio maps are within -10% 
of unity, indicating excellent agreement (see also Appendix of 
Konyves et al. 2010). 


4.2. Filamentary structure of the Aquila cloud complex 


As emphasized by [Men’shchikov et ah] ( |2010| ) and [Andre et ah] 
( 2010| ) and mentioned in Sect.[lJ filaments are widespread in the 
Herschel images of the Aquila region. Conceptually, an interstel¬ 
lar filament may be defined as any elongated structure in the ISM 
which is significantly denser than its surroundings. Eor the pur¬ 
poses of this paper, we adopt a minimum aspect ratio of - 3 and 
a minimum column density excess of - 10% with respect to the 
local background, i.e., >0.1, when averaged along 

the length of the structure. For more mathematical and algorith- 
mic de finit ions of a filament, th e reader is referred to jSousbiej 
( 2011| ) and [Men ’ shchikovj ( |20 13 j ) , respectively. 


In order to identify filaments in the high-resolution col¬ 
umn density map of the Aquila field, several methods were 
employed and compared. First, the contrast of elongated fea¬ 
tures was enhanced using a “morphological component analy¬ 
sis” (MCA) deco mposition of the m ap on a basis of curvelets 
and wavelets (e.g., |Starck et al.|2003| ). In such a decomposition, 
filamentary features are contained in the curvelet components, 
while roundish structures (e.g. dense cores) are contained in the 
wavelet components. Summing up all curvelet components led 
to the image shown in Fig. which provides a high-contrast 
view of the filaments after subtraction of core-like and other non- 
elongated structures (e.g. non-filamentary background). Given 
the typical filament width Wf^\ - 0.1 pc ( [Arzoumanian et ah] 
[2011[ ) and the relation Mune - So x Wfii between the central 
gas surface density So of a filament and its mass per unit length 
Mine (cf. Appendix A of [Andre et al.|2010| ), this curvelet com¬ 
ponent of the column density map is equivalent to a map of mass 
per unit length along the filaments. The white areas trace regions 
of the map where S x Wm is larger than half the critical value 
Mine,crit = (cf. |lnutsuka & Miyama||l997 ) and the fila¬ 
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ments are likely to be gravitationally unstable, i.e., supercritical^ 
with Miine ~ So X Wm > Miine,crit on the filament crest. 

A second, independent method used to trace filamentary 
structures in the mapped region was the multi-scale algo¬ 
rithm getfilaments ( [Men’shchikov[[2013[ ). Instead of tracing fil¬ 
aments directly in the observed images, getfilaments analyzes 
highly-filtered spatial decompositions of them (called “single¬ 
scale” images) across a wide range of scales (Sect. 2.3 of 
[Men’shchikov[[201^ . Using an automated iterati ve threshold¬ 
ing algorithm (Sect. 2.4.1 of Men’shchikov|2013[ ), getfilaments 
analyzes single-scale images and finds Icr intensity levels (on 
each spatial scale) that separate significant elongated structures 
from noise and background fluctuations. Setting to zero those 
pixels whose intensities are below the thresholds, the algo¬ 
rithm effectively “cleans” the single-scale images from noise and 
background. Fine spatial decomposition allows the algorithm 
to identify filaments as significantly elongated clusters of con¬ 
nected pixels on each spatial scale (Sect. 2.4.2 of [Men’shchikov[ 
[2013[ ), separating them from other (roundish) clusters of non- 
filamentary nature (e.g. sources or cores, noise peaks, isotropic 
backgrounds). Having produced the clean single-scale images of 
filamentary structures on each spatial scale, getfilaments recon¬ 
structs the intrinsic intensity distribution of the filamentary com¬ 
ponent of the images (largely free of sources, noise, and back¬ 
ground) by accumulating the clean decomposed images over all 
(or a range of) spatial scales (Sect. 2.4.3 of [Men’shchikov[2()T3 ). 
Finally, the algorithm generates mask images of filaments up 
to various transverse angular scales, as well as skeletons of the 
filament spines in th e format of fits images (see Sect. 2.4.4 of 
[Men’shchikov[[MT3| ). Filament extraction with getfilaments is 
fully automated and there are no free parameters involved. Fig¬ 
ure [^displays the filamentary network obtained by applying get¬ 
filaments to the high-resolution column density map shown in 
Fig.[^ For better visualization. Fig. shows a mask image cor¬ 
responding to elongated structures with transverse angular scales 
up to 320", equivalent to -0.4 pc at J = 260 pc. The color scale 
displayed within the filamentary mask corresponds to the col¬ 
umn density values in the input column density map (i.e., Fig.[^. 
The network of filaments outlined in this way (Fi g. [4| ) is very 
similar to that traced by the curvelet transform (Fig.[^. 

As a third, independent meth od to trace fila ments, we also 
applied the DisPerSE algorithnj^ ( Sousbie[201l| ). DisPerSE is a 
general tool to identify persistent topological features such as 
peaks, voids, and filamentary structures in astrophysical data 
sets. It traces filaments by connecting saddle points to max¬ 
ima with integral lines, following the gradient in a map. This 
method has already been used successfully to trace filamentary 
networks in Herschel images of nearby star-forming clouds (e.g.. 


Arzoumanian et al.[[2011l [Hill et al.[[2011 

Peretto et al.[[2012l 

Schneider et al.|2012[[Palmeirim et al.|2013 

). To trace filaments 


in the Aquila field, DisPerSE was run on the standard column 
density map (at 36.3"resolution) on a 6"/pixel scale where this 
pixel scale sets the resolution of the filament skeleton sampling. 
We used DisPerSE with a relative “persistence” threshold of 
4.8x10^^ cm“^, which corresponds to - 3 times the rms level 
of background column density fiuctuations in the low density 
portion (Ay ~ 2) of the column density image. “Persistence” is 
a measure of the robustness of topological features in the map 


^ Throughout this paper, by supercritical or subcritical filament, we 
mean a filament with a supercritical or subcritical mass per unit length 
(Miine > Mine,crit Or Miine < Mine,crit “ SCC ScCt. [^, rCSpeCtivcly. This 
notion should not be confused with the concept of a magnetically su¬ 
percritical or subcritical cloud/core (e.g.,[Mouschovias|1991[). 

^ see http://www2.iap.fr/users/sousbie/web/Jitml/mdexd41 d.html 
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Fig. 3. Left: Network of filaments in the Aquila cloud complex as traced by the curvelet transform compon ent (cf. Starck et'n!]|2003| ) of the 
Herschel high-resolution column density map shown in Fig.[^ Given the typical filament width Wfn ~ 0.1 pc ( [Arzoumanian et al.|2011| ) and the 
relation Mune ~ Eq X Wm between the centra l gas surface densit y Eq of a filament, this curvelet column density map provides information on the 
mass per unit length along the filaments (cf. [Andre et al.|2010j ), as indicated by the colo r bar on the right. The whit e areas highlight regions of 
the map where E x Wm exceeds half the critical mass per unit length Miine,crit = ^c^/G (cf. jinutsuka & Miyama|l997|) and the filaments are likely 
supercritical (Ep X Wfn > Mii ne,crit on th e fila ment crest). The overplotted blue skeleton marks the crests of the filaments selected with the DisPerSE 
algorithm of [Sousbi^ ( [2011[ ) (see Sect. [4.2[ for details). Right: Blow-up of the Aquila main subfield marked by the white square in the left panel, 
using the same color scale. 
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Fig. 4. Left: Mask of the filamentary network traced by getfilaments (m en’shchikov|2013| ) in the Herschel high-resolution column density map 
of the Aquila cloud complex. For better visualization, only angular scales up to 320" (i.e., ~0.4 pc at J = 260 pc) are shown. The color scale 
displayed wit hin the filament ary mask corresponds to colu mn de nsity values in the column density map (Fig.[^. The crests of the filaments traced 
by DisPerSE ( [Sousbie|2011[ ) are overlaid in blue (see Sect. [4.2[ for details). Right: Blow-up of the subfield marked by the white square in the left 
panel, using the same color scale. 


(see [Sousbie|2011| for details). Segments of filaments found by 
DisPerSE were assembled into longer filaments, with the con¬ 


straint that assembled segments did not form an angle larger than 
65°. The DisPerSE filaments were also trimmed to ensure that 
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the minimum column density along the resulting skeleton was 
5 X 10^^ cm“^ everywhere. This choice of DisPerSE parameters 
was adopted to facilitate the clean identification of dense, su¬ 
percritical filaments, which are most relevant to the problem of 
core formation and the present paper (see [Andre et al.|20TQ| and 
Sect. [5.3 1 below). From the output of DisPerSE, we then built a 
1-pixel-wide mask or skeleton image of the filament crests in the 
same way as [Arzoumanian et ar] ( |201 Ij ) did, after removing fila¬ 
mentary features shorter than 3 x 0.1 pc long (or ~ 80 pixels of 
3"). The resulting DisPerSE skeleton, which comprises a total of 
90 filaments, is overlaid in blue in both Fig.j^and Fig.|^ Owing 
to the adopted minimum column density, this DisPerSE skeleton 
is biased toward filaments which are either entirely or at least 
partly supercritical along their length. It nevertheless contains a 
dozen subcritical filaments. As can be seen by comparing Fig.[^ 
and Fig.[^ the three above-mentioned methods trace very simi¬ 
lar sets of filamentary structures. The agreement is particularly 
good in the case of supercritical filaments going over white areas 
in the image panels. 

The same filament profile analysis as described in [Arzouma^ 


nianet al.[ ( [2011[ ) was performed on the sample of filaments iden¬ 


tified here with DisPerSE, resulting in the distribution of filament 
inner widths shown in online Fig. A median FWHM width 
of 0.12 + 0.04 pc was found, which is very similar to the me¬ 
dian width reported by Arzoumanian et al. ( 2011[ ) for a smaller 
sample of 32 filaments in Aquila. 


4.3. Distribution of mass in the Aquila cloud 

Figure shows the probability density function (PDF) of col¬ 
umn density in the Aquila cloud complex as derived from the 
high-resolution column density map displayed in Fig. (Like¬ 
wise, the distribution of dust temperatures corresponding to the 
dust temperature map in Fig. is shown in online Fig. [T^.) 
The column density PDF is well fit by a log-normal distribu¬ 
tion at low column densities (i.e., 3 ^ Ay ^ 7) and by a power- 
law distribution at high column densities (i.e.. Ay ^7). Similar 
column density PDFs have already been reported in the liter¬ 
ature for other star-forming complexes based on near-infrared 
extinction data ( e.g.,[Kainulainen et al.|2009'] ) and Herschel ob- 
servations (e.g., [Schneider et al. [2013| ). As discussed by, e.g.. 


Mci = Mpixel/^H2 ^ 


extinction map of [Bontemps et al.[ ( [2010[ ) and with the total gas 
mass of -2.5x10"^ M© derived from the CO(l-O) observations of 
[Dame et ar| ( [2001[ ) over the same area. 

The same mass calculation was repeated for the pixels above 
a given column density, which led to the cumulative mass frac¬ 
tion of gas mass in the cloud as a function of column density 
shown in Fig.[^. For future reference, the fraction of dense gas 
mass above Ay =7-8 mag in visual extinction represents only 
-24-17% (5800-4200 M©) of the total cloud mass, respectively. 
A similar fraction of cl oud mass at Ay > 7 mag was reported by 
[Johnstone et al.[ ( [2004| ) in the case of the Ophiuchus main cloud. 
Clearly, the low (column) density regions in the map shown in 
Fig.[2account for most of the cloud mass. 

As the Aquila Rift lies quite close to the Galactic Plane, 
we have to consider the potential contamination of the column 
density map (Fig.[^ by background clouds along some lines of 
sight. To assess the importanc e of this potential c ontamination, 
we used the CO database of [Dame et al.[ ( [200T] ). Throughout 
the whole field shown in Fig. [IJ the most significant CO(l-O) 
emission was found in the same 5-7 kms“^ LSR velocity range. 
This correspondence suggests that the bulk of the CO emission 
comes from the same cloud complex at J - 260 pc. Two isolated 
patches with significantly higher LSR velocities (30-40 k m s“^), 
however, are present in the CO data of [Dame et al.[ ( [200r[ ), in the 
eastern corner and the relatively empty central part of the Her¬ 
schel field of Fig. respectively. Given that these patches are 
very local and contribute only modest CO emission, their infiu- 
ence on the above column density and mass estimates is minor. 
Based on the fraction of CO emission observed at high LSR ve¬ 
locities and the column densities derived from Herschel data, we 
estimate that the background patches cannot change the value of 
the total cloud mass given above - which excludes these patches 
- by more than 4%. We also stress that the exclusion of the CO 
high-VLSR areas from the distributions shown in Fig. has very 
little impact on the power-law slopes of the A-PDF and cumu¬ 
lative mass fraction plots since the background patches occupy 
only a small surface area. (Including the two patches in the N- 
PDF would change the power-law slope of -2.9 + 0.1 by much 
less than the quoted error bar.) The presence of these background 
clouds will nevertheless be taken into account when selecting 
dense cores belonging to the Aquila complex in Sect. [4.5[ 


Kainulainen et al. (2011), column density PDFs are a powerful 


tool to characterize molecular cloud structure and the transition 
from turbulence-dominated to collapsing, star-forming gas. The 
slope of the power-law distribution at high column densities can 
be readily related to the logarithmic slope of the equivalent ra¬ 
dial density profile expected in cloud collapse models (see, e.g., 
[Federrath & Kless^[2013[ ). The power-law tail of the PDF is 
particularly well developed in the case of the Aquila complex 
(see Fig.[^). The best power-law fit for Ay > 7 corresponds to 
JA/JlogAH2 oc 

The total mass of the Aquila cloud imaged with Herschel 
was derived from the column density map (Fig.[^ as: 


4.4. Multiwavelength core extraction with getsources 

Conceptually, a dense core is a single star-forming entity which 
may potentially form a star or a multiple system by gravitational 
collapse (e .g., Myers|19^[Ward-Thompson et al.|1994[ [Andre 


where Mpixei is the surface area subtended by one pixel at the 
adopted distance d = 260 pc of the cloud, fiu 2 = 2.8 is the mean 
molecular weight, mn is the hydrogen atom mass, and the col¬ 
umn density values in the map are summed up over all pixels. 
This procedure gave a total cloud mass of -2.4x10"^ M©. This es¬ 
timate of the total cloud mass from Aerschel data is in very good 
agreement with the total mass of -2.0x10"^ M© derived from the 


|et al.|2000} Di Francesco et al.|2007[ ). In practice, a core can be 

defined as the immediate vicinity of a local peak in the Herschel 
column density maps. In more mathematica l terms, a dens e core 
corresponds to a descending 2-manifold (cf. Sousbie[2011[ ) asso¬ 
ciated with a local peak in column density. This manifold defines 
a region in projection to the plane of sky whose map pixels are 
connected to the peak by lines following the gradient of the col¬ 
umn density distribution. 

To generate an extensive catalog of dense cores from RGBS 
data in the Aquila region, the parallel-mode SPIRE/PACS im¬ 
ages were processed with getsources, a multi-scale, multi- 
wavel ength source extraction algorithm ( [Men’shchikov et~ar 


2Q12| p This algorithm was designed primarily for extract¬ 
ing dense cores and young stellar objects (YSOs) in far- 


^ The RGBS first-generation catalog of cores presented in this pa¬ 
per (see Appendix A) was produced with the “November 2013” ma- 
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Fig. 5. a) Probability density function of column density (A-PDF) in the Aquila cloud, derived from the 18.2"-resolution column density image 
shown in Fig.^ The left and right axes give the actual and normalized numbers of independent beams per logarithmic bin in the column density 
map, respectively. (The right axis thus provides an estimate of a proper A-PDF whose integral is 1.) A powe r-law fit to the high column density tail 
at Ay > 7 gives dNld\ogNii 2 ^ The CO high-VlsR area in the eastern corner of the field (see Sect, 

b) Normalized cumulative mass fraction as a function of column density for the Aquila cloud (excluding t 


4.3 


has been excluded from this PDF. 


he CO high-VTsK area - see Sect. |4.3| ), 
based on the Herschel column density map shown in Fig.[^ The dense material to the right of the dashed vertical line (equivalent to Ay ~ 7-8) 
represents only ~17% of the total cloud mass, while the majority of the mass (~83%) corresponds to lower-density gas. A power-law fit to the 
cumulative mass fraction for Ay > 7 gives M{> oc ^ 


infrared/submillimeter surveys of Galactic molecular clouds 
with Herschel. The main features of the source extraction 
method, which may be conveniently divided into a detection 
and a measurement stage, can be summarized as follows (see 
[Men’shchikov et al.|2012[ for full details). 

At the detection stage, in contrast to the usual approach of 
detecting sources directly in the observed images, getsources 
analyzes ‘single-scale’ images (i.e., fine spatial decompositions 
of the original images - cf. Sect. |4.2[ ) across a wide range 
of scales and across all observed wavebands. This decompo¬ 
sition filters out irrelevant spatial scales and improves source 
detectability, especially in crowded regions and for extended 
sources. Usi ng an automated iterative thr esholding method (see 
Sect. 2.3 of |Men’shchikov et al.||201^ , getsources analyzes 
single-scale images and finds 3cr to 6cr intensity levels (on each 
spatial scale) that separate signals of significant sources from 
noise and background fluctuations. Setting to zero those pix¬ 
els whose intensities are below the thresholds, the algorithm ef¬ 
fectively ’cleans’ the single-scale images from noise and back¬ 
ground (including the filamentary component of the images). 
For detecting sources, getsources constructs a set of wavelength- 
independent single-scale detection images that preserve infor¬ 
mation in both spatial and wavelength dimensions (Sect. 2.4 of 
|Men’shchikov et ar]|2012| ). This multi-wavelength design com¬ 
bines data over all wavebands and thus naturally produces a 
wavelength-independent detection catalog with invariant source 
positions for all wavebands. Besides eliminating the need and 
problems of matching independent monochromatic extraction 
catalogs, the method also improves the detectability of weak 
sources and enables substantial super-resolution at wavelengths 
with lower spatial resolution. Sources are detected by getsources 
in the combined single-scale detection images by analyzing 
the evolution of their peak intensiti es and segmentation mask s 
across all spatial scales (Sect. 2.5 of |Men’shchikov et al.|2012| ). 


jor release of getsources (vl. 140127), which is publicly available from 
http://gouldbelt-herschel.cea.fr/getsources 


The spatial scale on which a source is brightest determines its 
characteristic size and corresponding footprint size. The latter is 
defined as the entire area that would give a non-negligible contri¬ 
bution to the integrated flux. The peak position of each source is 
determined from the wavelength-combined single-scale detec¬ 
tion images using the first moments of intensity (Appendix F 
of |Men’shchikov et al.||2012| ) measured over a range of spatial 
scales between the smallest scale on which the source appears 
and the characteristic scale on which the source is brightest. In 
effect, source coordinates are largely determined by the wave¬ 
bands with higher angular resolution and unaffected by large- 
scale emission. 

At the measurement stage, properties of detected sources are 
measured in the original observed images at each wavelength. 
These measurements go beyond simple aperture photometry 
since they are done together with background subtraction and 
deblending of overlapping sources (Sect. 2.6 of |Men’shchikov 
|et al.|[MT^ . Background is subtracted by linear interpolation 
under the source footprints found at the detection stage, con¬ 
strained by different angular resolutions in each waveband. The 
footprints must be at least as large as the beam size and their 
elongation must correspond to that of the source intensity distri¬ 
bution at that wavelength (Eq. 20 of |Men’shchikov et al.|2012| ). 
Overlapping sources are deblended in an iterative procedure 
that splits pixel intensity between blended sources, assuming 
a simple shape for their intensity distributions. The deblend¬ 
ing shape has a Gaussian-like circular pr ofile with somewhat 
stronger power-law wings (see Eq. 14 of [Men’shchikov et al. 
|2012| ) that should approximate the intensity profiles of observed 
sources. Local uncertainties of the peak intensities and integrated 
fluxes are given by the standard deviations estimated in elliptical 
annuli (covering areas of 20 observational beams) just outside 
the footprints. In crowded areas, the standard deviations are es¬ 
timated in expanded annuli outside of any of the overlapping 
sources (see Sect. 2.6 of |Men’shchikov et al. 112012] ). Aperture 
corrections are applied by getsources using tables of the encir- 
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cled energy fraction values for the actual PSFs provided by the 
PACS and SPIRE ICCs ( |Balog et al.|20T4l|Bendo et al.|2013| ). 

Source extraction with getsources is fully automated and 
there are no free parameters involved: default configuration pa¬ 
rameters have been extensively tested and fine-tuned to work in 
most practical cases. For the production of the ‘first-generation’ 
catalogs of starless and protostellar cores from the RGBS, the 
following two-pronged extraction strategy has been adopted. 
Two sets of dedicated getsources extractions are performed, op¬ 
timized for the detection of dense cores and YSOs/protostars, 
respectively. 

In the first set, all of the Herschel data tracing column density 
are combined at the detection stage, to improve the detectability 
of dense cores. The detection image is thus combined from the 
clean 160yum, 250yum, 350yum, and 500yum maps, together with 
the high-resolution column density image (see Sect. |4.1| ) used as 
an additional “wavelength”. The latter is added to the combined 
detection image to ensure that detected sources correspond to 
genuine column density peaks. Furthermore, the 160 yum compo¬ 
nent to the detection image is “temperature-corrected” to reduce 
the effects of strong, anisotropic temperature gradients present 
in parts of the observed fields, such as in the vicinity of the W40 
HII region in Aquilaj^ The temperature-corrected 160 yum map is 
obtained by converting the original observed 160yum map (13.5" 
resolution) to an approximate column density image, using the 
color-temperature map derived from the intensity ratio between 
160yum and 250yum (at the 18.2" resolution of the 250yum map). 
Simulations on synthetic emission maps including model cores 
(see, e.g.. Sect. |4 8 [ below and Appendix jB.lj ) confirm the valid¬ 
ity of this approach to detecting dense cores. 

A second set of getsources extractions is performed to trace 
the presence of self-luminous YSOs/protostars and discriminate 
between protostellar and starless cores. Here, the only Herschel 
data used at the detection stage come from the 70yum image. 
Indeed, the presence of point-like 70 yum emissio n traces the in- 
temal luminosity of a protostar very well (e.g., [Dunham et al.j 
2008| ), and HerscheljVACS observations of nearby (J < 500 pc) 
clouds even have the sensitivity to detect candidate “first hydro¬ 
static cores”, the very first and lowest-luminosity stage of proto¬ 
stars (cf. [Pezzuto et al.|2()T^ . 

At the measurement stage of both sets of extractions, source 
properties are measured at the detected positions of either cores 
or YSOs/protostars, using the observed, background-subtracted, 
and deblended images at all five Herschel wavelengths, plus the 
high-resolution column density map. The advantage of this two¬ 
pronged extraction strategy is that it provides more reliable de¬ 
tections and measurements of column-density cores and 70yum 
luminous YSOs/protostars, respectively. 

4 . 5 . Selection and classification of reliable core detections 

Here, we summarize the criteria adopted to select various types 
of dense cores from the raw source lists produced by the two 
sets of multi-wavelength getsources extractions described at the 
end of Sect. [4.4[ For each source type, the following prescribed 
criteria should be met at the same time. 

^ In the presence of an anisotropic radiation field, due to a closeby 
HII region for instance, radiative transfer calculations show that the far- 
infrared emission expected from a starless core at, e.g., 160 pm is not 
centered on the column density peak but is shifted toward the source 
of illumination. Using a “temperature-corrected” IbOyum map instead 
of the original 160yum map at the detection stage in getsources allevi¬ 
ates this problem and helps to better trace the intrinsic position of the 
underlying column density core. 
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Selection of candidate dense cores (either starless or protostel¬ 
lar) from the ''core” set of extractions 

- Column density detection significance greater than 5, where 
detection significance here refers to a single-scale analog 
to a classical signal-to-noise ratio (S/N) [see Eq. (17) of 
Men’shchikov et al. 2012] in the high-resolution column 
density map; 

- Global detection significance over all wavelengths [see 
Eq. (18) of Men’shchikov et al. 2012] greater than 10; 

- Global “goodness” > 1, where goodness is an output qual¬ 
ity parameter of getsources, combining global signal-to- 
noise ratio and source reliability, and defined in Eq. (19) of 
Men’shchikov et al. (2012); 

- Column density measurement signal-to-noise raticj^ (S/N) 
greater than 1 in the high-resolution column density map; 

- Monochromatic detection significance greater than 5 in at 
least two bands between 160 pm and 500 pm\ and 

- Flux measurement with S/N > 1 in at least one band between 
160 pm and 500 pm for which the monochromatic detection 
significance is simultaneously greater than 5. 

Selection of candidate YSOs from the "protostellar” set of 
extractions 

- Monochromatic detection significance greater than 5 in the 
70 pm band; 

- Positive peak and integrated fiux density at 70 yum; 

- Global “goodness” greater than or equal to 1; 

- Flux measurement with S/N > 1.5 in the 70 pm band; 

- FWHM source size at 70 pm smaller than 1.5 times the 
70 pm beam size (i.e., < 1.5 x 8.4" or < 12.6"); and 

- Estimated source elongation < 1.30 at 70 pm, where source 
elongation is defined as the ratio of the major and minor 
FWHM sizes. 

The discussion of the Herschel-i&eniified sample of protostars 
and YSOs in Aquila will be presented in a complementary paper 
(Konyves et al., in prep.; see Maury et al. 2011 for a preliminary 
subsample around W40 and Serpens-South). 

Selection of candidate starless cores and protostellar cores 

- After cross-matching the selected dense cores with the can¬ 
didate YSOs/protostars, a selected dense core is classified as 
’starless’ if there is no candidate 70 pm YSO within its half¬ 
power (high-resolution) column density contour. 

- Conversely, a selected dense core is classified as ’protostel¬ 
lar’ if there is a candidate 70 pm YSO within its half-power 
column density contour. 

- The most reliable SED of a selected protostellar core is ob¬ 
tained by combining the 70 pm fiux density from the “pro¬ 
tostellar” extractions with the 160 pm, 250 pm, 350 pm, and 
500 pm fiux densities from the “core” extractions. 

Post-selection checks 

All of the cores automatically selected according to the above 

The measurement S/N is esti mated at the measurement step of the 
getsources extractions (see Sect. \4A\ and characterizes the flux mea¬ 
surement uncertainties. In crowded situations, the measurement S/N of 
a source with a high “detection significance” at the detection step can 
be low because of large deblending and background-subtraction uncer¬ 
tainties. 
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criteria were visually inspected in the SPIRE/PACS and column 
density images (see blow-up maps in Figs. |A.3| & |A.4| ). Any du¬ 
bious source was removed from the final catalog of cores pre¬ 
sented online Table A.l (see below). 

To eliminate from our discussion of Aquila cores extragalac- 
tic contaminants that may be misidentified as cores or YSOs, we 
also cross-matched all selected sources with the NASA Extra- 
galactic Databas 4^ (NED), but no close match (within 6") was 
found. 

Likewise, we checked likely associations between the se¬ 
lected Herschel cores and objects in the SIMBAD database or 
the combined c2d and Gould Belt Spitzer database (Dunham et 
al. 2013; Allen et al., in prep.). Any matches are reported in the 
online catalog (Table A.l). In particular, 27 associations with a 
Spitzer source were found using a 6" matching radius. 

In the eastern corner of the field shown in Fig. there are 
two known dense clumps (ISOSS J18364-0221 SMM1/SMM2) 
with > 30 km/s LSR velocities from molecular line measure¬ 
ments ( [Birkmann et al.|20'Q6| ). Their LSR velocities correspond 
to a kinematical distance of ~2.2 kpc. We therefore excluded 
from our Aquila discussion 23 candidate prestellar cores and 6 
protostellar cores (shown as yellow triangles in Fig. [T]) lying in 
the hig h-VLS R CO area of the Herschel field mentioned at the end 
of Sect. |4.3| These cores are nevertheless listed (with appropriate 
comments) in the online catalogs. 

In the post-selection phase, we also used another source 
extraction method to generate an “alternative-algorithm” fiag 
for our getsources master source catalog entries. For this pur¬ 
pose, we used CSAR (Cardiff Sourcefinding AlgoRithm - |Kirk,| 
|J. M. et al.|2013| ), a hierarchical source-finding algorithm, which 
we applied to the high-resolution column density map. In the 
Aquila entire field, the fraction of matches between selected 
starless cores and CSAR-detected sources is ~45%, adopting 
a matching separation of less than 6" between peak positions. 
The resulting ’CSAR’-fiag appears in the online catalogs to in¬ 
dicate if a given getsources core was independently detected 
by CSAR. Several reasons explain the relatively low fraction of 
matches with CSAR-detected sources based on peak positions. 
First, CSAR is a very conservative source-finding algorithm, es¬ 
pecially in crowded areas where the multi-scale nature of get¬ 
sources makes it easier to detect blended objects. Second, CSAR 
is a ‘monochromatic’ algorithm which detects sources in a sin¬ 
gle observed image (here the column density map, intrinsically 
noisier than the Herschel images) and does not benefit from the 
multi-wavelength design of getsources (significantly improving 
the detectability of weak sources). Third, the present core sam¬ 
ple is dominated by starless cores which tend to have fiat-topped 
density profiles ( |Ward-Thompson et al.||1994'| ) and whose peak 
positions are not as well defined as the peak positions of YSOs. 
Relaxing our matching condition somewhat, we note that ~ 65% 
of all selected cores include the peak position of a CSAR source 
within their FWHM contours. Moreover, as much as ~ 85% of 
the surface area occupied by the FWHM ellipses of our selected 
cores in the high-resolution column density map is also covered 
by the outer footprints of CSAR sources. Based on our visual in¬ 
spection checks, we are confident that all getsources detections 
reported in online Table A.l are robust. (The reader can judge by 
looking at the blow-u p ma ps prov ided online for each source - 
see examples in Figs. |A.3| & |A.4[ ). For the sake of robustness at 
this ‘first-generation’ stage, ~ 20% of the automatically-selected 
getsources cores were visually rejected for reasons such as map 
artifacts, sources seen only at some Herschel bands but not in 
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the column density image, or sources seen only in the column 
density map but not at any of the Herschel bands. 

Our getsources selection and classification procedure re¬ 
sulted in a final sample of 709 Aquila cores (not counting 40 ob¬ 
jects - 34 starless cores, including 23 candidate prestellar cores, 
and 6 protostellar cores - in the high-VLSR CO area), comprising 
651 starless cores and 58 protostellar cores. A total of 446 star¬ 
less cores were classified as candidate prestellar cores on the ba¬ 
sis of their locations in a mass versus size diagram (see Sect. |4.7| 
and Fig. [7] below). The spatial distribution of the prestellar and 
protostellar cores is shown in Fig. overplotted on the high- 
resolution column density map of the cloud. 

The observed properties of all selected cores are given in 
the accompanying online catalog (cf. Table A.l in Appendix A). 
The contents of Table A.l are as follows: core running number 
(Col. 1), RGBS source name (Col. 2), J2000 equatorial coor¬ 
dinates (Cols. 3 & 4), then, for each Herschel wavelength, de¬ 
tection significance (5), peak fiux density and error [(6) & (7)], 
contrast over local background (8), peak fiux density in a 36.3" 
beam (9), total integrated flux density and error [(10) and (11)], 
major & minor FWHM diameters [(12) & (13)], and position an¬ 
gle of the major axis (14), followed by detection significance in 
the high-resolution column density map (15), peak H 2 column 
density at 18.2" resolution (16), column density contrast over 
local background (17), peak column density in a 36.3" beam 
(18), column density of local background (19), major & minor 
FWHM diameters and position angle of the major axis in the 
high-resolution column density map [( 20 ), ( 21 ), & ( 22 )], number 
of Herschel bands at which the core is significant (23), ’CSAR’- 
flag (24), core type (25), SIMBAD counterpart(s) (26), Spitzer- 
c2d counterpart if any (27), and comments (28). 


4.6. Derived core properties 


The SED fitting procedure used to derive core properties was 
similar to the procedure described in Sect. |4.1| for the produc¬ 
tion of the column density map. Here, the SEDs were con¬ 
structed from the integrated flux densities measured by get¬ 
sources for each extracted core (see Fig. |A.1| ) and the SED data 
points were weighted by l/cr^^, where (Terr corresponds to the 
flux measurement error estimated by getsources for each point. 
(In contrast to Sect. |4.1| where the dominant source of error 
was the calibration uncertainty, the errors on source flux es¬ 
timates are primarily driven by uncertain background subtrac¬ 
tion.) The modified blackbody fits to the ob served SEDs were 
performed with the MPCURVEFIT routine ( |Markwardt||2009| ) 
in IDF. These SED fits provided direct estimates of the mass 
and line-of-sight-averaged (SED) dust temperature for most of 
the selected cores. The core masses were derived assuming the 
same dust opacity law as in Sect. |4.1| and a distance d = 260 pc 
for the Aquila complex. The angular FWHM size estimate re¬ 
turned by getsources for each core (as measured at 18.2" res¬ 
olution in the high-resolution column density map) was con¬ 
verted to a physical core radius assuming the same distance. 
Two estimates of the core radius are provided (see online Ta¬ 
ble A.2). The first estimate is a deconvolved radius, calculated as 

^deconv = (FWHmI^^ where FWHM and HPBW 

denote the physical sizes corresponding to the FWHM angular 
size of the core and the HPBW resolution of the high-resolution 
column density map, respectively. The second estimate simply 
corresponds to the observed average FWHM size of the core (ge¬ 
ometrical average between the major and minor FWHM sizes). 
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Fig. 6. Distributions of beam-averaged column densities (left) and beam-averaged volume densities (right) at the resolution of the SPIRE 500 pm 
observations for the population of 446 candidate prestellar cores in Aquila (solid curves). In both panels, the dashed curves show the corresponding 
distributions for all 651 selected starless cores. 


In principle the first value provides a more accurate estimate of 
the intrinsic core radius, but it is affected by significantly larger 
uncertainties than the second value in the case of marginally re¬ 
solved cores. In the case of a self-gravitating prestellar core, both 
values provide estimates of the core outer radius under the as¬ 
sumption that such a core can be appro ximately desc ribed as a 
critical Bonnor-Ebert (BE) sphere (e.g., |Bonnor|I956| ). (Indeed, 
a critical BE sphere of outer radius /?be has a column density 
profile approaching that of a Gaussian distribution of EWHM 
diameter ~ /?be-) A peak (or central beam) column density, an 
average column density, a central-beam volume density, and an 
average volume density were also derived for each core based on 
its estimated mass and radius. The central-beam column density 
was estimated from the peak flux densities of the core at the res¬ 
olution of the SPIRE 500 pm observations (HPBW = 36.3" or 
~ 0.046 pc at J = 260 pc) using an SED fitting procedure similar 
to that described in Sect. |4.I| The central-beam volume density 
no (at the same resolution) was derived from the central-beam 
column density No assuming a Gaussian spherical distribution, 
for which no = A^o/( V^cr), where cr is the standard deviation 
of the Gaussian distribution. The distributions of column densi¬ 
ties and volume densities for the population of starless cores are 
shown in Eig.[^ 

All of the derived properties are provided in online Table A.2 
for the whole sample of selected Herschel cores. The contents of 
Table A.2 are as follows: core running number (Col. 1), RGBS 
core name (Col. 2), J2000 equatorial coordinates (Cols. 3 & 
4), deconvolved and observed core radius (Cols. 5 & 6), esti¬ 
mated core mass and corresponding error (Cols. 7 and 8), SED 
dust temperature and corresponding error (Cols. 9 and 10), peak 
column density at 36.3" resolution (Col. 11), average column 
density measured before and after deconvolution (Cols. 12 and 
13), beam-averaged peak volume density at 36.3" resolution 
(Col. 14), average volume density derived before and after de- 
convolutio n (C ols. 15 and 16), Bonnor-Ebert mass ratio (Col. 17 
- see Sect. |4.7| ), core type (Col. 18), and comments (Col. 19). 

Since color correction factors are usually small, we did not 
apply any color corrections to the measured fiux densities. In¬ 
stead, like [Kelly et ^ ( |20I2| ), we adopted an elevated calibra¬ 
tion uncertainty representing multiple sources of uncertainties. 
Our adopted calibration uncertainties for the SED data points 
were 10-20% for the PACS 70-160 pm bands and 10% for the 
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SPIRE 250/350/500 pm bands, respectively, which are conser¬ 
vative values compared to the HSC-recommended point source 
calibration uncertainties^ 

The robustness of the SED fits was assessed by using and 
comparing two successive runs of the fitting routine with slightly 
different weighting schemes for each source. In the first run the 
70 pm data point was included in the fit and the getsources de¬ 
tection errors were used to weigh the SED data points, while 
in the second run the 70 pm point was not fitted and the (more 
conservative) measurement errors were used to weigh the SED 
data points. The detection errors on significant data points were 
typically ~ 15% (comparable to the absolute calibration un¬ 
certainty), while the median measurement errors ranged from 
~ 30% to ~ 70% depending on wavelength (being typically 
higher at I60yum). The results of the SED fits were accepted for 
a given source if I) significant fiux measurements exist for this 
source in at least three Herschel bands, 2) the source has a larger 
integrated fiux density at 350 pm than at 500 pm, and 3) there 
was less than a factor of 2 difference between the core mass es¬ 
timates derived from the two fit runs. About 68% of the starless 
cores had reliable SED fits. The corresponding distribution of 
SED dust temperatures is shown in online Eig. [T^. Compari¬ 
son with the distribution of dust temperatures in the background 
cloud (online Eig. [T^) indicates that the Aquila starless cores 
are somewhat colder than the parent cloud, as expected (cf. [Royj 
jet al.|20I4D . 

The masses of the starless cores for which the SED fit results 
were rejected were directly estimated from the measured inte¬ 
grated fiux density at the longest significant wavelength in each 
case, assuming optically thin dust emission at the median dust 
temperature found for starless cores with reliable SED fits (i.e., 
11.5 + 2 K outside the W40/Sh62 areas and 14.5 + 3 K within the 
higher radiation field areas W40 and Sh62). The corresponding 
cores have more uncertain properties and are marked as having 
“unreliable SED fits” in the last column of online Table A.2. 

Accuracy of the core mass estimates 

Uncertainties in the dust opacity law alone induce uncertainties 


The photometric point-source calibration uncertainty is less than 7% 
for the PACS bands (Balog et al. 2014) and ~ 5% for the SPIRE bands 
(Bendo et al. 2013). 
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in low-m a ss (and intermediat e-mass) starless cores (e.g., |My 
1983[ I Andre et al.||2007| ). For each object, we estimated 


ers 


the thermal BE mass, Mbe, from the deconvolved core radius 
^deconv measured in the high-resolution column density map (see 


of up to a fac tor ~ 1.5-2 in the core mass estimates. As men¬ 
tioned in Sect. |4.1| the dust opacity law adopted here and in other 
RGBS papers, namely = 0.1 x (A/300jdm)~^ cm^Ig, is likely 

appropriate to better than 50% in the 160-500 fim range for col- _ 

umn densities between ~ 3 x 10^^ cm“^ and ~ 10^^ cm“^ (cf. Sect. |4.6| ) assuming a typical gas temperature of 10 K. 
|Roy et al.|2014| ). 

In addition to the dust opacity, another systematic effect af¬ 
fects the accuracy of our simple SED mass estimates. A single¬ 
temperature graybody fit to the integrated flux densities can 
only provide an average value of the dust temperature for each 
source and neglects any variation in the dust temperature within 
the source. In reality, starless dense cores, which are externally 
heated objects, are known from both radiative transfer calcula¬ 
tions (e.g., [Evans et al.|2001t |Stamatellos et al.|2007|) and, e.g., 

Herschel observations (e.g., [Nielbock et al. 201 2} [Roy et al.j 
2014| ) to have a stratified temperature structure with a signifi¬ 


cant drop in dust temperature toward core center. In such a sit¬ 
uation, the average dust temperature derived from a global SED 
fit can sometimes significantly overestimate the mass-averaged 
dust temperature within a starless core, leading to an underes¬ 
timate of the core mass. The magnitude of this effect is very 
modest (< 20%) for low column density cores such as B 68 ( jRoyj 
jet al.||2014| ) but increases to up to a factor ~ 2 for high-density 
cores with average column densities ^ 10^^ cm“^. In the case 
of spatially-resolved cores with good signal-to-noise data, tech¬ 
niques such as the Abel-inversion method ( |Roy et aL]|2014| ) or 
the COREEIT method ( [Marsh et~ar]|2014| ) can help to retrieve 
the intrinsic temperature structure and derive more accurate mass 
estimates. We did not attempt to use such techniques here. Based 
on the results of the simulations performed to estimate the com¬ 
pleteness of the survey (see Sect. [4.8[ below and Appendix [B.l[), 
however, we estimate that the SED masses listed in Table [A.2l 
for starless cores are likely underestimated by ~ 20-30% on av¬ 
erage compared to the intrinsic core masses, mainly due to the 
fact that the SED dust temperatures tend to slightly overestimate 
the intrinsic mass-averaged temperatures of starless cores. The 
column densities and volume densities listed in Table IA.2I and 
used in Eig.|^(see also Eig.j^below) have not been corrected for 
this small effect. 


4.7. Selecting self-gravitating prestellar cores 

Conceptually, a dense core is deemed to be prestellar if it is 
both starless and sel f-gravitating (cf. [Andre et al.[[200'0} [Di 


[Erancesco et al.[[2007t[Ward-Thompson et al.[[2007 ). Such star¬ 

less cores will most likely form (proto)stars in the future. Lack¬ 
ing spectroscopic observations for most of the Herschel cores, 
we used the thermal value of the critical Bonnor-Ebert (BE) mass 
(4^BE,crit - Bonnor 1956) to assess whether a starless core was 
self-gravitating or not based on the value of its BE mass ratio 
(Tbe = 4^BE,crit/4^obs- The Critical BE mass can be expressed as 


4^be,c] 


2.4 Rbe cJG, 


where Rbe is the BE radius, Cg the isothermal sound speed, and 
G the gravitational constant. In the presence of significant non- 
thermal motions, the thermal BE mass should be replaced by 
a modified BE mass obtained by substituting the total (ther¬ 
mal -h non thermal) one-dimensional velocity dispersion for the 
isothermal sound speed in the above formula. The simplified ap¬ 
proach adopted here, where the nonthermal component of the 
velocity dispersion is neglected, is justified by observations of 
nearby cores in dense gas molecular tracers such as NH 3 and 
N 2 H^ lines, which show that nonthermal motions are negligible 



Fig. 7. Mass versus size diagram for the entire population of 651 starless 
cores identified with Herschel in the Aquila cloud. The core FWHM 
sizes were measured with getsources in the high-resolution column 
density map (Fig.[^ and deconvolved from an 18.2" (HPBW) Gaus¬ 
sian beam; the vertical dashed line marks the corresponding physical 
HPBW resolution at 4 = 260 pc. The core masses were derived as ex¬ 
plained in Sect. [4.6[ Typical error bars are shown for both well-resolved 
and marginally-resolved cores (to the right and the left of the vertical 
dashed line, respectively). The 292 robust prestellar cores (for which 
Q^be < 2 - see text) are shown as filled blue triangles, the other (candi¬ 
date) prestellar cores as open blue triangles, and the rest of the starless 
cores as open green triangles. The red curve shows the empirical lower 
envelope used to select the 446 candidate prestellar cores in the dia¬ 
gram (i.e., QfBE < 5 X {HPBW^^^IFWHM^^^)^"^ - see text), based on 
the Monte-Carlo simulations described in Sect. [4.8[ For comparison, 
models of critical isothermal Bonnor-Ebert spheres at T = 7 K and 
r = 20 K are plotted as black solid lines. The mass-size corr elation 
observed for diffuse CO clumps ( [Elmegreen & Falgarone|199^ is dis¬ 
played as a shaded yellow band. The blue curve marks a column density 
level corresponding to 5 x A^H 2 ,rms. where A^H 2 ,rms is the typical rms level 
of cirrus noise fluctuations at Ay ~ 7 in the column density map (see 
Fig. B.l in Appendix B). The arrow at the upper right indicates how 
the Herschel cores (and the blue curve marking the cirrus noise level) 
would move in the diagram using a distance of 415 pc instead of 260 pc 
(see Appendix [C]). 

In practice, we used the positions of the Herschel cores in a 
mass versus size diagram (Fig. [7]) to distinguish between candi¬ 
date prestellar cores and unbound starless cores, after deriving 
a reasonable lower envelope for sel f-gravitating cores in such 
a diagram. In our first-look papers ( [Konyves et al.[[20Tm [An- 
[dre et al.| 2010 [ ), the criterion adopted to define this lower enve¬ 
lope was simply ctbe < 2 , by analogy with the usual criterion 
to select self-gravitating objects based on the virial mass ratio 
(ofyir = Mvir/Mobs < 2 - e.g., [Bertoldi & McKee[1992] ). Adopting 
the same criterion here led to a first sample of 292 robust prestel¬ 
lar cores, shown as filled blue triangles in Fig. [7| However, the 
Monte-Carlo simulations performed in Sect. [4.8[ toow to assess 
the completeness of the survey suggest that this criterion may be 
too restrictive, in the sense that it selects only ~ 85% of the simu¬ 
lated BE cores detected by getsources after source classification. 
Eor this reason, we also derived a less restrictive lower enve- 
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lope (shown as a red curve in Fig. [7]) based on the results of our 
Monte-Carlo simulations. This second, empirical lower envelope 
contains > 95% of the simulated BE cores after getsources ex¬ 
traction, and corresponds to the following, size-dependent limit¬ 
ing BE mass ratio: a^E < 5 x (HPBW^^^jFWHM^^ where 
FWHM^^^ is the measured FWHM source diameter m the high- 
resolution column density map and = 18.2" is the 

HPBW resolution of the map. The limiting BE mass ratio varies 
from ~ 2 for well-resolved cores with FWHM^^^ ~ 0.1 pc to 


~ 5 for unresolved cores with FWHM^^ ~ HPBW ^^^. The 
reason why one has to be more flexible and use a larger limiting 
BE mass ratio for unresolved or marginally resolved cores is that 
the intrinsic core radius (and therefore the intrinsic BE mass) is 
more uncertain for such cores. 

Based on the latter criterion, 446 of the 651 starless cores 
in the Aquila entire field were classified as candidate prestellar 
cores. All of the 292 robust prestellar cores belong to the wider 
sample of 446 candidate prestellar cores. These two samples of 
cores refiect the uncertainties in the classification of detected 
starless cores as gravitationally bound or unbound objects, which 
are fairly large for marginally-resolved cores. In the absence of 
higher-resolution observations, the status of the 155 candidate 
prestellar cores which do not match the first criterion (ckbe ^ 2 ) 
is more uncertain (these cores a re ma rked as “tentative bound” in 
the last column of online Table |A.2] ). We will thus consider both 
samples of prestellar cores in the discussion presented in Sect. 
below. 

The mass vs. size distribution of the entire population of 
selected starless cores (Eig. |7]) shows a spread of deconvolved 
FWHM sizes between ~ 0.01 pc and ~ 0.1 pc and a range in 
core mass between ~ 0.03 M© and ~ 10 M©. The high frac¬ 
tion of self-gravitating cores (~ 45% or ~ 69%, depending on 
whether the robust or the candidate sample is adopted) is re¬ 
flected in the locations of the Aquila starless cores in this mass 
vs. size diagram. The selected robust prestellar cores are clus¬ 
tered around (or above) the mass-size relations expected for crit¬ 
ical BE isothermal spheres with gas temperatures between 7 K 
and 20 K (parallel black solid lines in Eig.|7]). Besides, they are 
more than an order of magnitude denser than typical CO clumps 
(yell ow band in Eig. [7|), which are mostly unbound structures 
(e.g., [Elmegreen & Ealgarone|1996[[Kramer et al.|1998| ). 


4.8. Completeness of the prestellar core survey 

To estimate the completeness of our census of prestellar cores 
in Aquila, we performed Monte-Carlo simulations (see Ap¬ 
pendix |BT]). We first constructed clean maps of the background 
emission at all Herschel wavelengths (including a column den¬ 
sity plane), by subtracting the emission of the compact sources 
identified with getsources. We then inserted a population of 
~ 5600 model Bonnor-Ebert-likj^ cores throughout the clean- 
background images to generate a full set of synthetic Herschel 
and column density images of the region. The model cores were 
given a fiat mass distribution (dA/dlogM oc from 0.02 M© 

to ~ 30 M© and were assumed to follow sl M oc R mass versus 
size relation appropriate for isothermal spheres. The dust con¬ 
tinuum emission from the synthetic Bonnor-Ebert cores in all 



0.1 1.0 10.0 100.0 
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Fig. 8. Completeness curve of our Herschel sample of candidate 
prestellar cores as a function of true core mass (solid line), as estimated 
from the Monte-Carlo simulations described in Sect. |4.8| For compari¬ 
son, the dashed line shows the global completeness curve predicted by 
the model discussed in Appendix |B.2| 


Herschel bands was simulated using an extensive grid of spher¬ 
ical dust radiative tra nsfer models constructed by us with th e 
MODUST code (e.g., [Bouwman et~ST 2000^ |Bouwman||2QQl| ). 
Compact source extraction for several sets of such synthetic 
skies was performed with getsources in the same way as for the 
observed images. 


Based on the results of these simulations (see Appendix |B.1| ) 
for further details), we estimate that our Herschel census of can¬ 
didate prestellar cores is > 90% complete above a true core 
mass of -0.3 M© which corresponds to an observed core mass 
of -0.2 Mq on average, given that observed masses are typi¬ 
cally underestimated by - 20-30% due to the internal tempera- 

and Fig. |B.2| 

in Appendix]^. Likewise, our sample of robust prestellar cores 
is estimated to be - 80% complete above a true core mass of 
-0.3 Mq or an observed core mass of -0.2 M©. The complete¬ 
ness curve of the Aquila core survey as a function of true core 
mass is plotted in Eig.[^ 


ture structure of starless cores (see end of Sect. 4.6 


In reality, the completeness level of the core survey is ex¬ 
pected to be background dependent. In an effort to assess the 
magnitude of this dependence, we constructed a simple model 
of the prestellar core population and core extraction process de¬ 
scribed in Appendix |B.2| This model shows that the complete¬ 
ness of prestellar core extractions does decrease as backgr ound 
cloud column density and cirrus noise increase (see Eig. |B.6| ) 
but suggests that the global completeness curve of the prestellar 
core survey in Aquila is consistent with that inferred from our 
Monte-Carlo simulations (compare the dashed and the solid line 
in Eig.[^. 


We use the term “Bonnor-Ebert-like” because the model cores were 
given the density structure of critical isothermal Bonnor-Ebert spheres, 
but their dust temperature distributions resulted from radiative transfer 
calculations and were thus not strictly isothermal, in agree ment with 
detail ed observational studies of individual cores (see, e.g., [Roy et al.| 


2014 for the example of B68). 


Armed with a good understanding of the completeness of the 
core survey, we discuss in Sect.j^below the global properties of 
the dense core population and their connection with the filamen¬ 
tary structure of the cloud complex on the basis of statistically 
representative observational results. 
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5. Discussion 


5.1. Lifetimes of Herschel prestellar cores 

As our Herschel survey provides an essentially complete cen¬ 
sus of prestellar cores in the Aquila cloud, the core statistics can 
be used to set constraints on the typical lifetime of prestellar 
cores and the timescale of the core formation process. Following 
a technique introduced by |Beichman et al] ( 1986| ) in the con¬ 
text of IRAS data, a rough estimate of the lifetime of prestellar 
cores can be obtained by comparing the number of starless cores 
found with Herschel to the number of Class II YSOs detected by 
Spitzer in the same region. The underlying assumptions are 1) 
that all starless cores will evolve into YSOs in the future, and 2) 
that star formation proceeds at a roughly constant rate, at least 
when averaged over an entire cloud. In the ~ 3° x3° field covered 
by Herschel, and excluding the dubious, small area with higher 
LSR velocities in the eastern corner of the column density map 
(see Fig. and Sect. |4.3| ), our survey revealed a total of 651 
starless cores, including 446 candidate and 292 robust prestellar 
cores, while the combined c2d and Gould Belt S pitzer surveys 
detected 622 Class II YSOs ( [Dunham et al.|2013| Allen et al., in 
prep.). Adopting a reference lifetime of 2 Myr for Class II YSOs 
( Evans et~aL]|2009| ), these numbers lead to typical lifetimes of 
~ 2 Myr, -1.4 Myr, and ~ 0.9 Myr for the global populations of 
Herschel starless core^ candidate prestellar cores, and robust 
prestellar cores, respectively. Given the large sizes of the popu¬ 
lations of starless cores and YSOs in Aquila, the main sources of 
error in these estimates come from the fact that some starless or 
even candidate prestellar cores may never evolve into YSOs, as 
they may be “failed cores” that will disperse before collapsing 
(e.g. jVazquez-Semadeni et al.||2005] ), and from the uncertainty 
in the number and lifetime of Class II YSOs in Aquila. Com¬ 
bining the constraints coming from the two samples of observed 
prestellar cores, our best estimate of the global lifetime of the 
prestellar core phase is ^pre = 1.2 + 0.3 Myr. 

We have a large enough sample of cores in Aquila to inves¬ 
tigate a possible trend between core lifetime and core density. 
Figure shows a plot of estimated core lifeti me versus average 
volume density, similar to that introduced by jJessop & Ward-j 
[Thompson] ( [2000| ), but for the sample of Herschel-i&cn\ihc& 
candidate prestellar cores in Aquila. In this plot, the Aquila data 
are represented by blue triangles and compared to literature data 
(black crosses) from [Ward-Thompson et al.[ p007[ ). The blue 
solid line and filled triangles represent the estimated trend be¬ 
tween core lifetime and average core density, where the latter 
quantity was obtained by dividing the observed mass of each 
core by the deconvolved estimate of its volume [i.e., core den¬ 
sity reported in col. (16) of online Table A.2]. As can be seen, 
the plot suggests that the typical lifetime of prestellar cores de¬ 
creases from -1.4 Myr for cores with average volume density 
^10"^ cm“^ to a few times 10"^ yr for cores with average volume 
density ^ 10^-10^ cm“^. Moreover, the estimated core lifetimes 
lie between one free-fall time (%, lower dashed line in Fig.]^, the 
timescale expected in free-fall collapse, and lOx^ff (upper dashed 
line in Fig. [^, roughly the timescale expected f or highly sub- 
critica l cores undergoing ambipolar diffusion (e.g., [Mouschovias[ 
1991 ). At the median average volume density - 4 X 10"^ cm of 


the candidate prestellar cores identified with Herschel, the es¬ 
timated core lifetime is - 0.75 Myr or -4-5%. The densest 


The lifetime estimate quoted for starless cores is given under the as¬ 
sumption that presently unbound starless cores are still growing in mass 
and will become gravitationally bound and prestellar in the futu re (cf. 
Simpson et al.|2011 Belloche et al. 2011 Kirk, J. M. et al.|20T^ . 



Fig. 9. Plot of estimated lifetime against minimum average volume den¬ 
sity (blue solid line and filled triangles) for the population of 446 candi¬ 
date prestellar cores identified with Herschel in the Aquila cloud (blue 
triangles), similar to the “JWT” plot introduced by [Jessop & Ward-[ 


Thompson (2000). The error bars only reflect sfN counting uncertain¬ 
ties. Literature data from [Ward-Thompson et al.| ( [2007[ ) are shown as 
black crosses for comparison. The two parallel dashed lines correspond 
to the free-fall timescale (%) and a rough approximation of the ambipo¬ 
lar diffusion timescale (10 x %). 



10 ‘ 1 10 * 
Derived core mass, M [Mq] 


Fig. 10. Plot of average volume density versus observed core mass for 
the sample of 446 candidate prestellar cores. The blue solid triangles 
mark the median deconvolved volume density for each mass bin. (For 
comparison, the black open triangles show the median 36.3"-beam- 
averaged densities for the upper three mass bins.) The error bars corre¬ 
spond to the interquartile range of densities in each mass bin. The data 
points become very uncertain at the high-mass end due to the small 
number of cores in the higher mass bins (see core mass function in 
Fig. below). (No interquartile range can be plotted for the last bin 
which contains only two cores.) Note the weak, but significant, correla¬ 
tion between core density and core mass above - 2-3 Mq. The red curve 
represents a parabolic fit to the data points. 


cores in our sample, which have beam-averaged volume den¬ 
sities ^ 2 X 10^ cm“^ at the resolution of the 500 pm data and 
average deconvolved volume densities ^10^ cm“^, have a much 
shorter lifetime - 0.02-0.05 Myr or - %, suggesting they may 
evolve essentially on a free-fall timescale. Indeed, the tentative 
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Background column density* [10^^ cm"^] Background column density, Nh^ [lO^^cm 


Fig. 11. Left: Distribution of background cloud column densities for the population of 446 candidate prestellar cores identified with Herschel in 
the whole Aquila field. The vertical dashed line marks a fiducial threshold at an equivalent visual extinction level A^ ~ 7 mag (cf. Heiderman 


|et al.||2010[|Lada et al.|2010| | Andre et al.||2010||2014| ). Right: Normalized “probability” of finding a Herschel prestellar core as a function ot 
background column density (blue histogram with error bars), obtained by dividing the number of prestellar cores detected with Herschel below a 
given background column density level by the total surface area covered by the RGBS survey below the same level. The red curve shows a simple 
fit of the form PcoreC^v) = 1 - exp(a xA^ + b) with a - -0.17 and b - 0.86. The vertical dashed line marks Ay ~ 7 mag as in the left panel. 


presence of a power-law tail in the distribution of beam-averaged 
core densities above ~ 10^ cm“^ (see Fig.[^) suggests that these 
cores may be undergoing nearly free-fall collapse. 

In this context, it is worth pointing out that density may 
not be the only relevant parameter and that core evolution may 


also be mass dependent as suggested by, e.g., Hatchell & Fuller 
( 2008| ). Indeed, we observe a weak positive correlation between 
core density and core mass above ~ 2-3 Mq (see Fig. [TO, indi- 


.bg 


and is such that most (-90% prestellar cores are 
found above a background column density corresponding to 
Ay^ - 7 and a background gas surface density - 150 Mo/pc^. 
As already emphasized by [Andre et al.| ( |MT0||2014| ), the shape 
of the distribution shown in Fig. [11^ strongly supports the exis¬ 
tence of a column density threshold for the formation of prestel- 


More precisely, 88% of the candidate prestellar cores and 92% of the 


robust prestellar cores lie at A^ > 1 
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lar cores. The existence of such a threshold had been suspected 
for a long time, based on the results of ground-based millimeter 
and submillimeter surveys for cores in, e.g., the Taurus, Ophi- 
uchus, and Perseus clouds (e.g., jOnishi et al.||1998t [Johnstone 


eating that the most massive prestellar cores in our sample tend 
to be the densest objects. Assuming that the lifetime of a core 
is proportional to its free-fall time, the correlation in Fig. [T^ 
suggests that prestellar cores more massive than - 2-3 Mq may 
evolve on significantly shorter timescales than the majority of 
the cores in our sample, which have masses ~0.1-2 Mq (see 
Fig. below). This finding would be consistent with the re¬ 
sults of earlier searches for high-mass prestellar cores (i.e., pre¬ 
cursors to stars > 8 Mq) in massive star-forming regions which 
have shown that such cores, if they exist at all, are extremely 
rare with lifetimes comparable to (or shorter than) the free-fall 
timescale ( [Motte et al.|2007| ). 

5.2. Evidence of a column density threshold for prestellar 
core formation 

Figure shows the distribution of background cloud col¬ 
umn densities entire population of 446 candidate 

prestel lar c or es id entified with Herschel in the Aquila cloud (see 
Sects. [4.4[ & [4.7[). This distribution shows a steep rise above 


et al.[[2004t[Kirk, H. et al.[[2006[ ). These early claims, however, 

were not completely convincing due to the limited column den¬ 
sity sensitivity and spatial dynamic range of ground-based obser¬ 
vations, hence their limited capability to probe prestellar cores 
and the parent background cloud simultaneously. The Herschel 
results presented in this paper provide a much stronger case 
for a (column) density threshold. We stress that the distribution 
of cl oud mass as a function of column density (see Fig. in 
Sect. [4.3[ ) and the background-dependent completeness level of 
our survey for prestellar cores make the threshold even more sig¬ 
nificant than Fig. [TT^ suggests. Indeed, ~ 85% of the mass in 
the Aquila cloud is at column densities lower than Ay ~ 7 (see 
Fig. [^) and ~ 95% of the surface area covered by the Her¬ 
schel survey is below Ay ~ 7. Furthermore, the completeness 
level of our Herschel census for prestellar cores is not limited 
by sensitivity (as was typically the case for earlier ground-based 
surveys), but by “cirrus confusion noise” (see Appendix B), and 
is better in Ay <7 areas than in Ay >7 areas (see Fig. [B.6| ). 
Therefore, if prestellar cores were distributed randomly in the 
cloud, we would be much more likely to detect prestellar cores 
in Ay <7 areas than in higher column density regions. Fig¬ 
ure [TT^ already shows that this is clearly not the case. To further 
strengthen the point, we plot in Fig.[TT]) a probability function of 
finding a prestellar core as a function of background column den¬ 
sity, obtained by normalizing the number of prestellar cores de¬ 
tected below a given background column density by the total sur¬ 
face area imaged by Herschel below the same background col¬ 
umn density level (for a related probability function in the case 
of t he submm continuuni cores detected by SCUBA in Perseus, 
Hatchell et al. 2005). The probability function, P®Q^^g(Ay), 


see 


shown m Fig. 115 increases by more than an order of magni¬ 
tude between Ay - 4 and Ay - 10, and looks like a smooth 
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step function. It is very well fit by the simple exponential step 
function /^core(^v) = 1 “ exp(-0.17 x Ay + 0.86). 


5.2.1. Comparison with models of the star formation rate 

There is some debate in the literature as to whether the kind of 
results shown in Fig. [TT] reflect a true column density threshold 
for star formation or whether the efficiency of the star forma¬ 
tion process simply increases gradually with (column) density 
(cf. Hatchell et al. |2005| ). Starting with the work of |Krumholz| 
|& McKee] ( |2005| ), a number of theoretical models of the star 
formation rate (SFR) in molecular clouds have been proposed 
based on the general idea that star formation is regulated by in¬ 
terstellar turbulence and that clouds typically convert 6ff ~ 1% 
of their molecular gas mass into stars per (local) free-fall time 


(e.g.,|Padoan & Nordlund|[20TH Hennebelle & Chabrier||2011 


IKrumholz et al. 120121- see also 


'ederrath & Klessenl20i2land 


Padoan et al.||2014| for overviews and comparisons of the mod¬ 
els). In the “multi-freefall” versions of these theoretical m odels. 


which are most appropriate to fit real observationj^(cf., Hen 


[nebelle & Chabrier]|2Qllt [Federrath & Klessen||2Q12| ), there is 
not necessarily any sharp (column) density threshold, but the 
SFR drops significantly at low densities because of a signifi¬ 
cant increase in the local free-fall time (see also the related dis¬ 
cussion by |B urkert & Hartmann|20 1 3] ) . In Fig.[T^ we compare 
the observed core formation efficiency (CFE) as a function of 
background column density with the prediction of the simplified 
multi-freefall model of [Hennebelle & Chabrier|(2011| ). Here, we 
define the observed core formation efficiency as CFEobs(Av) = 
AMpores(AV)/ AM h on h(AV) whcrc AMcores(^v) is the mass of the 
prestellar corej^identified with Herschel in a given bin of back¬ 
ground Ay values and AMcioud(^v) is the cloud mass estimated 
from the Herschel column density map in the same Ay bin. In 
the multi-freefall model, the fraction of gas mass converted into 
core mass per unit time is simply — x —where 6ff ~ 1% (see 

^core 

above), 6core ~ 40% is the star formation efficiency at the level 
of an individual prestellar core (see Sect. |5.5| below), and ts(p) 
is the local free-fall time at the local gas density p. Over the typ¬ 
ical lifetime of prestellar cores ^pre ~ 1 Myr (see Sect. 1^, the 
expected core formation efficiency is thus: 


CFEmff (p) = - ^ w \ • 

^core %(P) 

In order to use this formula, we had to estimate the local gas den¬ 
sity and free-fall time in the Aquila cloud. To do so, we made 
use of the fact that the cloud surface area above a given col¬ 
umn density level S(> N 112 ) scales as the column density PDF 
shown in Fig.[^ and in particular features a well-defined power- 
law tail S (> A^H 2 ) ^ ^H 2 ^ high column densities (Ay > 5-7). 
In spherical geometry, this is indicativq^ of a power-law den- 


In the initial, “single-freefall” model of 


Krumholz & McKee 


(20051 


the relevant timescale is the free-fall time evaluated at the mean d ensity 
of the cloud, tff(po), and there is no density dependence at all. Hen-j 
[nebelle & Chabrier| ( |2011| ) and [Federrath & Klessen[ ( |201^ have 5 ' 


that this model generally underestimates the SFRs determined by [Hei-[ 
derman et al. ( 2010j ) in nearby clouds. 


Both AMcoresfAy) and AMcioud(Av) represent observed masses di¬ 
rectly estimated from the Herschel data using SED dust temperatures 
and the dust opacity of assumptions given in Sect.|^ AMcores(Av) was 
not corrected for the small ~ 25% effect due to the fact that the SED 
mass values tend to slightly underest imat e the intrinsic core masses ac¬ 
cording to our simulations (see Sect.[4.6[and Appendix [B.lj l. 

In the case of a cloud with a spherical radial density distribution, p oc 
it is easy to show that both the column density PDE, dN/d\ogNii 2 ^ 



Fig. 12. Plot of the observed differential core formation efficiency (CPE) 
as a function of background column density expressed in Ay units (blue 
histogram with error bars), obtained by dividing the mass in the form 
of prestellar cores in a given column density bin by the cloud mass 
observed in the same column density bin. The right axis gives the cor¬ 
responding star formation rate per unit gas mass (in units of yr"^), esti¬ 
mated assuming a local star formation efficiency 6core = 0-4 at the core 
level and a prestellar core lifetime tpre = 1.2 Myr (see text). Por compar¬ 
ison, the black curve shows the prediction of the multi-freefall version 
of the turbu lence regulated model of the star formation rate initially 
proposed by [Krumholz & McK^ ( |2005j ) (see [Hennebelle & Chabrier[ 
[2011| ). The vertical dashed line marks the same fiducial threshold at 
Ay ~ 7 as in Pig.f^ The horizontal dotted line marks the rough asymp¬ 
totic value of the CfiE ~ 15% at Ay > 15, corresponding to a star for¬ 
mation rate per unit gas mass ~ 5 x 10"^ yr"k 


sity distribution p oc r~^-^ for the dense gas and is consistent 
with large-scale cloud contraction above Ay > 5-7. Under the 
assumption of a roughly spheroidal ambient cloud, we then de¬ 
rived the effective volume density, and effective free- 

fall time, ^ff(Ay), of the gas as a function of background cloud 
density expressed in Ay units. Applying the above multi-freefall 
formula, this allowed us to obtain the core formation efficiency, 
CFEniff(Ay) = ^ X predicted by the multi-freefall model 
as a function of Ay, for direct comparison with CFEobs(Ay). 
As can be seen in Fig. the Herschel observations indicate 
a much sharper transition than the multi-freefall model does, 
between a regime of negligible prestellar core formation effi¬ 
ciency at Ay < 5 and a regime of roughly constant CFE -15% 
at Ay > 15. Furthermore, we stress that differential comp lete- 
ness between low and high column density areas (see Fig. [B.6[ 
in Appendix B) implies that the real transition between the two 
regimes is in fact somewhat sharper than indicated by the blue 
histogram in Fig.[^ On this basis, we argue for the presence of 
a true physical threshold for prestellar core formation around a 
fiducial value Ay ^1, although the observed transition is clearly 
not infinitely sharp like a true Heaviside step function. 

Interestingly, a very similar extinction threshold at A^^^^ - 7 
has independently been observed with Spitzer in the spatial dis- 
tribution of YSQs in nearby c louds ( [Heiderman et al.[2010[[Lada 
[et al.|2010|[Evans et al. 2014[- see also Sect. [5.6[ below). Follow¬ 
ing [Andre et al.[ ( 2010 [2014| ), we interpret this star formation 
threshold in terms of the quasi-universal filamentary structure of 
molecular clouds in Sect. l5.4l below. 


and the surface area, S (> N uj), scale as AH 2 " 
Federrath & Klessen 2013| >. 


where m = (see, e.g.. 
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Fig. 13. Comparison of the spatial distribution of the prestellar core population identified in Sect. 4.7 using getsources (blue triangles) with the 
footprints of all the filaments traced with DisPerSE (red contours), both overlaid on the curvelet component of the high-resolution column density 
map (cf. Fig. [^. The gray scale corresponds to the color scale of Fig. The red contours outline 0.1 pc-wide footprints around the crests of 
filaments. The green contours correspond to Ay = 7 in the column density map smoothed to a resolution of 5'. 


5.3. Spatial distribution of Herschel cores and connection 
with filaments 


As already pointed out in earlier RGBS papers (e.g., Andre et al. 
|2010[ [Men’shchikov et aL]|201Q| ), there is a very close corre¬ 
spondence between the spatial distribution of compact dense 
cores and the network of filaments identified in the Herschel 
column density map of the Aquila cloud. Furthermore, candi¬ 
date prestellar cores and embedded protostars are preferentially 
found within the densest filaments with supercritical masse s per 
unit length (i.e., Mune > Miine,crit = 2 c^/G - see Sects. [^& |4.2[ ) 
(e.g., [Andre et al.|2010[|2014| ). 


The connection between cores and filaments is illustrated in 
Figs.[^& l^and can be quantified in detail based on the cen¬ 
sus of cores presented in Sects. |4.4| |4.5| & [4.7[ and the census 
of filaments described in Sect. f4.2[ To this end, a mask image 
of the filament “footprints” was constructed by convolving the 
filamentary skeleton traced with DisPerSE and shown in Fig. 
and Fig. [^ with a Gaussian kernel corresponding to a typical fil¬ 
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ament inner width ~ 0.1 pc ( [Arzoumanian et aL]|2011| ), i.e., an 
angular width ~ 80" at the distance of Aquila. An alternative 
mask image of the filaments, similar to that shown in Fig.[^ was 
created by considering all transverse angular scales up to 80" 
in the multi-scale decomposition performed by getfilaments (see 
Sect. [4.2[ ). The core positions were then compared with these two 
sets of 0.1-pc filament footprints to estimate the fraction of cores 
associated with filaments. The results of this comparison, sum¬ 
marized in Tableindicate that a very high fraction 
of prestellar cores are closely associated with filaments, i.e., lie 
within 0.1-pc filament footprints. This correspondence is illus¬ 
trated in Fig. where the 0.1-pc footprints of the filaments 
traced with DisPerSE (cf. Fig. [^ are outlined by red contours 
and the population of 446 candidate prestellar cores identified 
with getsources in Sect. [^ are superimposed as blue triangles. It 
can be seen that most (~ 70%) of the candidate prestellar cores 
lie within the red filament footprints. Likewise, ~ 80% of the 
robust prestellar cores lie within the red filament footprints. A 
more detailed view of this connection in the Aquila “main sub- 
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Fig. 14. Upper: Curvelet component (cf. Starck et al. 2003 1 ) of a portion 
of the Herschel high-resolution column de nsity map shown in Fig. \T\ 
Given the typical filament width of ~0.1 pc ( [Arzoumanian et al.|20li| ), 
this map is equivalent to a map of the mass per unit length along the fil¬ 
aments (cf. [Andre et al.|2010] ), as indicated by the color bar on the right. 
The white areas highlight regions of the map where the filaments have a 
ma ss per unit length larger tha n half the critical value Miine,crit = 2cl/G 
(cf. [Inutsuka & Miyama|1997[ ) and are thus likely to be gravitationally 
unstable (see Sect. [4.2[ and Fig.[^. The contours overlaid in blue out¬ 
line the 0.1 pc-wide footprints of the filaments traced with DisPerSE in 
Sect. [4.2[ (cf. red contours in Fig.[^. Lower: Same map as in the upper 
panel with the locations of candidate prestellar cores and protostellar 
cores overlaid as blue triangles and green triangles, respectively. 


field”, including Serpens-South and W40, is provided by Fig. 14 


Table 1. Fractions of cores associated with filaments in Aquila. 



DisPerSE 

getfilaments 


0.1 pc 0.2 pc 

0.1 pc 0.2 pc 

All filaments: 

prestellar ON-fil. 
starless ON-fil. 

71%-78% 

60% 

81%-88% 

72% 

83%-87% 

75% 

84%-89% 

77% 


Supercritical segments: 


prestellar ON-fil. 

66%-75% 

76%-84% 

76%-81% 

77%-83% 

starless ON-fil. 

55% 

66% 

67% 

69% 


Notes. The upper part of this table gives the fractions of prestel- 
lar/starless cores found inside the 0.1 pc and 0.2 pc-wide filament foot¬ 
prints constructed with DisPerSE and getfilaments over the Aquila en¬ 
tire field (see text). The lower part of the table provides similar cores 
fractions when only supercritical portions of the filaments are consid¬ 
ered. Here, for the sake of simplicity, a portion of a filament was clas¬ 
sified as either supercritical or subcritical based on whether the local 
column density in the clean background column density image (after 
subtracting the contribution of cores with getsources) was equivalent to 
^d,back ^ 7 Qj. ^d,back ^ 7 resp ectively, assum ing a constant filament 
width ~ 0.1 pc (see Sect. [4.2[ and Sect. \5A\ . The lower fractions of 
the ranges quoted for prestellar cores correspond to candidate prestellar 
cores, the higher fractions to robust prestellar cores. 


which shows the locations of candidate prestellar cores overlaid 
on the curvelet component (cf. [Starck et~ar][20Q3 1 of the high- 


resolution column density map (see Andre et al. 


[2010[ for an 


early version of the same view). It is important to stress that the 
connection between cores and filaments does not strongly de¬ 
pend on the precise definition adopted for a filament or on the 
algorithm used to trace filaments. In particular, as can be seen 
in Table [TJ the values found for the fractions of cores associated 
with filaments using getfilaments footprints are very similar to 
the values found using DisPerSE footprints. 

Table also reports the fractions of cores found within su¬ 
percritical portions of filaments. For the sake of simplicity, in 
the present paper focusing primarily on cores, our classifica¬ 
tion of filament segments as either supercritical or subcritical 
relies on the assumption of a constant filament width ~ 0.1 pc 
( [Arzoumanian et ar][2011[ ) and is based on the local column 
density measured in the clean background column density im¬ 
age (after subtracting the contribution of cores with getsources). 
To take into account the fact that the transverse column den¬ 
sity profiles of supercritical filaments featu re power-law wings 
which extend beyond the 0.1 pc inner width ( [Arzoumanian et al" 


2011[ [Palmeirim et al.[[2013[ ), we also considered 0.2 pc-wide 


filament footprints and provide corresponding core fractions in 
Table For example, the well-studied Serpens South filament 
and Taurus-B211/B213 filament have transverse column density 
profiles which extend up to ~0.4-0.5 pc in radius on average and 
equivalent widths of ~ 0.2 pc ( [Hill et al.|2012t [Palmeirim et al 


[2013[ ). While a detailed discussion of the radial column density 
profiles of the present filament sample is out of the scope of this 
paper, simple comparison of the line masses obtained by integra¬ 
tion over the filament profiles with the line masses derived from 
the central column densities of the filaments using a characteris¬ 
tic inner width of 0.1 pc suggests that the equivalent widtlf^of 
the supercritical filaments traced here with DisPerSE (see end of 
Sect. [4.2[ ) is also typically ~ 0.2 pc. 


Here, we define the equivalent width of a filament as the effective 
width Weff such that the line mass integrated over the filament profile is 
= 2^0 X Ikeff, where So is the central surface density of the filament. 
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5.4. Mass budget in the cloud and interpretation of the star 
formation threshold in terms of the filamentary structure 

Our Herschel census of prestellar cores and filaments allows us 
to derive a detailed mass budget in the Aquila cloud. Below the 
fiducial column density threshold at ~ 7, ~ 10%-20% of 
the gas mass is in the form of (mostly subcritical) filaments and 
< 1% of the cloud mass is in the form of prestellar cores. Above 
^back _ 7 _ 50%-60% of the cloud mass is in the form of 
(mostly supercritical) filaments and a fraction /pj-e ~ 15% + 5% 
of the mass is in the form of prestellar cores. We note that /pre 
roughly corresponds to the asymptotic core formation efficiency 


value reached at > 15 in Fig. 12 The fraction of cloud mass 
in the form of filaments reaches a very high value ~ 75% above 
^back _ 10. In attempt to quantify further the relative contribu¬ 
tions of cores and filaments to the cloud material as a function 
of column density, we compare in Fig. the column density 
PDFs observed for the cloud before any component subtraction 
(blue histogram, identical to the PDF shown in Fig.[^, after sub¬ 
traction of dense cores (red solid line), and after subtraction of 
both dense cores and filaments (black solid line). To generate this 
plot, we used getsources to create a column density map of the 
cloud after subtracting the contribution of all compact cores, and 
gethlaments to construct another column density map after also 
subtracting the contribution of filaments. Although there are ad¬ 
mittedly rather large uncertainties involved in this two-step sub¬ 
traction process, the result clearly suggests that filaments domi¬ 
nate the mass budget of the Aquila cloud at high densities. 



Fig. 15. Comparison of the global column density PDF in the Aquila 
cloud (blue histogram with statistical error bars, identical to that shown 
in Fig.[^ to the column density PDF measured after subtraction of dense 
cores (red solid line) and the PDF measured after subtraction of both 
dense cores and filaments (black solid line with statistical error bars). 
The vertical dashed line marks the same fiducial threshold at Ay ~ 
7 as in Fig. [TT] This plot illustrates that filaments make up a dominant 
(~ 50%-75%) fraction of the dense gas mass at Ay > 7-10, and that 
dense cores contribute only a small ( ^ 15%) fraction of the dense gas 
(except perhaps at the very highest column densities). 


Since filaments appear to make up a dominant fraction of the 
dense gas material at > 7 within which the vast majority of 
prestellar cores are observed (see Fig. [D, and since the spatial 
distribution of prest ellar cores is strongly correlated with fila¬ 
ments (see Sect. |5.3| ), it is tempting to interpret the star formation 
threshold discussed in Sect. |5.2| in terms of the quasi-universal 
filamentary structure of molecular clouds (cf. [Andre et al.|2014] ). 
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Given the typical width Wm ~ 0.1 pc measured for filaments 
(Arzoumanian et al. 2011) and the relation Mune ~ x be¬ 
tween the central gas surface density Eq and the mass per unit 
length Miine of a filament (cf. Appendix A of [Andre et al.|2010| ), 
the threshold at A^^^^ ~ 7 or -150 M© pc“^ corresponds 
to within a factoip^ of < 2 to the critical mass per unit length 
Afiine,crit = 2 cJ/G - 16 M© pc“^ of nearly isothermal, long cylin¬ 
ders (see Inutsuka & Miyama 1997) for a typical gas temperature 
T - 10 K. Thus, the prestellar core formation threshold approxi¬ 
mately corresponds to the threshold above which interstellar fil¬ 
aments become gravitationally unstable ( [Andre et al.|2010[ ). 


5.5. Prestellar CMF and link with the IMF 


The prestellar core mass function (CMF) derived from the sam¬ 
ples of 446 candidate and 292 robust prestellar cores identified 
in the whole Aquila cloud (see Sect. [4.7[), excluding the CO 


& [4.5[ ), is shown in the form of a differential mass distribution 
in Fig. (see dark blue histograms and light blue shade). The 
mass distributi on o f the wider sample of 651 starless cores se¬ 
lected in Sect. [4.5[ is plotted as a green histogram for compar¬ 
ison. The 90% completeness level of our Herschel census of 
preste llar c ores, as estimated both from Monte-Carlo sim ulatio ns 
(Sect. [4.8[ ) and the simple model described in Appendix [B.2[ is 
marked by the vertical dashed line. We stress that the differen- 
tial C MF presented here (see [Konyves et al.[[2010[ [Andre et^ 
[2010[ for preliminary versions of this CMF) is based on a core 
sample ~ 2-9 times larger than the CMFs derived from earlier 
ground-based studies (e.g.,|Motte et al.[[ 1998 [ [Johnstone et al. 
[2000t[Stanke et al.][2006t[Mves et al. [[2007[[Enoch et al.[[2008[ ) 

and that its shape is therefore much more robustly defined. In 
particular, it suffers very little from the arbitrary choice of mass 
bins, a well-known disadvantage of differential mass functions 
(e.g., [Reid & Wilson[2006[ ), except perhaps at the very high mass 
end (e.g. at M ^ 5 M©, where the number of cores per mass bin 
drops to less than 10 in Fig. [T^. Note also that, while we pre¬ 
ferred to display the differential form of the CMF in Fig. [T^ be¬ 
cause it is more intuitive and easier to compare with the IMF, we 
used the cumulative form - which is independent of binning and 
thus amenable to cleaner statistical tests - to quantify the resem¬ 
blance of the observed CMF to several well-known functional 
forms. 


high-VLSR area in the eastern comer of the field (see Sects. 4.3 


As can be seen in Fig.[^ the Aquila prestellar CMF is well 
fit by a lognormal distribution (solid and dashed red curves for 
the samples of robust and candidate prestellar cores, respec¬ 
tively) and very similar in shape to the system IMF advocated 
by [ChabrTer| ( [2005) . Performing a non-parametric Kolmogorov- 
Smirnov (K-S) test (see, e.g. [Press et ar][1992[ ) on the corre¬ 
sponding cumulative mass distributions N(>M) indicates that the 
observed prestellar CMF is statistically indistinguishable at the 
97% confidence level from a lognormal mass function with cen¬ 
tral mass 0.45 + 0.2M© and standard deviation 0.52 + 0.05 above 
the completeness mass limit ~ 0.2 M©. For comparison, the log¬ 
normal part of the [Chabrier[ ( [2005[ ) system IMF has a central 
mass of 0.25 M© and a standard deviation of 0.55 in logioM. The 
error on the two parameters of the lognormal fit to the prestel- 


Strictly speaking, the formal agreement between x Wfn and 
Mine,crit IS cven better than 10%. For several reasons (e.g., factor of - 2 
spread in filament width and distribution of filament inclination angles), 
however, the column density threshold is not a sharp boundary bu t a 
smooth transition (see discussion in Sect. 6.2 of [Andre et al.|2014| ), as 
also observed in Fig. 
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Fig. 16. Differential core mass function (dA/dlogM) of the 651 star¬ 
less cores (dark green histogram), 446 candidate prestellar cores (up¬ 
per blue histogram and open triangles), and 292 robust prestellar cores 
(lower blue histogram and filled triangles) identified with Herschel in 
the whole Aquila field. The error bars correspond to ^lN statistical 
uncertainties. The shaded area in light blue reflects the uncertainties 
in the prestellar CMF arising from the uncertain classification of ob¬ 
served starless cores as gravitationally bound or unbound objects (see 
Sect. |4.7| ). The dashed blue histogram and open squares show how 
the prestellar CMF would change at the high-mass end after correc¬ 
tion for a possible differential timescale bias (see text). The 90% com¬ 
pleteness level of the prestellar core sample is indicated by the verti¬ 
cal dashed line (see Sect. |4.8| and Appendix B). Lognormal fits to the 
CMF of robust and candidate prestellar cores (solid and dashed red 
curves, respectively), as well as a power-law fit to the high-mass end 
of the CMF (black solid line) are superimposed. The two lognormal fits 
peak at 0.62 M© and 0.34 M©, and have standard deviations of ~ 0.47 
and ~ 0.57 in logioM, respectively. The power-law fit has a slope of 
-1.33 + 0.06 (compared to a Salpeter slope of -1.35 in this format) . 
The IMF of single stars (correct ed for binaries - e.g., |Kroupa||2001| ), 
the IMF of multiple systems (e. g.,|Chabrier|2005]), a nd the typical mass 
distribution of CO clumps (e.g., Kjamer et al.|1998| ) are also shown for 
comparison. 


lar CMF (i.e., central mass and standard deviation) are mainly 
driven by the uncertain classification of observed starless cores 
as gravitationally bound or unbound objects, which leads to two 
slightly different CMF shapes for the samples of robust and can¬ 
didate prestellar cores (blue shaded area in Fig. [^. The high- 
mass end of the Aquila CMF above I M© is also consistent with 
a power-law mass function, dA^/dlogM oc - 33 + 0 . 06 ^ ^ 

significance level of 87%. Here, the error bar on the power-law 
exponent was derived from the range of values for which the K- 
S significance level is larger than 68% (corresponding to I cr in 
Gaussian statistics). This function is very similar to the Salpeter 
power-law IMF which is dA^/dlogM oc in this format. (We 

note, however, that given the limited range of core masses probed 
by our data a power law does not provide a significantly bet¬ 
ter fit to the high-mass end of the CMF than a pure lognormal 
fit.) In contrast, the CMF observed above I Mq differs from the 
shallower power-law mass distribution of CO clumps and cl ouds 
(dA^/dlogM oc - e.g. Blitz||l993 Kramer et al.||l998 ) at a 


very high confidence level. The probability that the CMF can be 
consistent with dA^/dlogM oc M-O-’ is only Pk-s ~ 7.7 X lO”’. 

A possible caveat to the similarity between the Salpeter IMF 
and the prestellar CMF at the high-mass end should be men¬ 
tioned, however. As pointed out by [Clark et al.| ( |2007] ), if cores of 


different mass evolve on different timescales then the observed 
CMF may not be representative of the intrinsic prestellar core 
mass function. This is because an observer is more likely to de¬ 
tect long-lived cores than short-lived cores. Therefore, if there is 
a correlation between core lifetime and core mass, then the ob¬ 
served CMF can be significantly distorted compared to the “ini¬ 
tial” prestellar core mass distribution. In the present sample of 
prestellar cores, there is essentially no correlation between core 
density and core mass below ~ 2 Mq but a weak positive corre¬ 
lation above ~ 2-3 Mq (Fig.[^, suggesting that prestellar cores 
more massive than ~ 2-3 Mq may evolve to protostars somewhat 
faster than lower mass cores do (see end of Sect. |53- To quan¬ 
tify the importance of this potential differential timescale bias 
on the CMF, we have overplotted in Fig. [T^ a weighted ver¬ 
sion of the observed CMF (blue open squares and dashed his¬ 
togram), obtained by weighting the number of prestellar cores 
observed in each mass bin by a factor inversely proportional to a 
mass-dependent free-fall time. The latter was estimated for each 
mass bin by using the parabolic fit to the observed correlation 
between core density and core mass shown by the red curve in 
Fig. [T^ As can be seen in Fig. [T^ the weighted CMF is in¬ 
distinguishable from the unweighted CMF for M < 2Mq, but 
somewhat shallower above ~ 2Mq. (A K-S analysis indicates 
that the high-mass end of the weighted CMF above 2 Mq is con¬ 
sistent with a power-law mass function, dA^/dlogM cx 
at a K-S significance level of 90%.) The main effect of the dif¬ 
ferential timescale correction is to broaden the prestellar CMF, 
leaving the peak mass at ~ 0.5 Mq essentially unchanged. 

As already discuss ed by [Andre et al. ( [20I0[ ) and [Konyves 
[et al.[ ( [20T0| ) (see also [Alves et al.[[2007 ), the observed CMF 
is consistent with an essentially one-to-one mapping between 
prestellar core mass and stellar system mass^ i.e., M*sys = 
6 core X ^core). where 6core represents the efficiency of the con¬ 
version process from core mass to stellar system mass, i.e., the 
star formation efficiency within an individual prestellar core. The 
peak of the prestellar CMF is at 0.45 + 0.2 Mq in observed core 
mass, suggesting a real peak at 0.6 + 0.2 Mq in terms of intrinsic 
prestellar core mass, after correcting the observed masses up¬ 
ward by ~ 25% due to the fact that the SED mass values tend to 
slightly underestimate the intrinsic core masses according to our 
simulations (see Sect. [4. 6 [a nd Appendix [B.I[ ). Our data therefore 
suggest that 6core ~ 0-4^ • 

It is also interesting to investigate possible variations in the 
CMF as a function of local cloud environment, in particular de¬ 
pending on whether the cores lie within or outside dense fila¬ 
ments. In LI64I (Orion A), for instance, [Poly chroni et al.[ ( [20T3] ) 
reported that the cores lying on filaments were generally more 
massive than those lying off filaments. Figure [T^ compares the 
CMF derived for the candidate prestellar cores lying on filaments 
(fight blue histogram) to the CMF of the candidate prestellar 
cores lying off filaments (magenta histogram) and to the global 
prestellar CMF in Aquila (upper dark blue histogram). It can be 
seen that the prestellar CMF observed on filaments is very simi¬ 
lar to the global prestellar CMF. A two-sample K-S test confirms 
that these two CMFs are indistinguishable at a > 95% confidence 
level. On the other hand, there is a marginal indication that the 


As pointed out by a number of authors (e.g. Delgado-Donate et al. 


[2003 1 [Goodwin et al.[2008[[Hatchell & Fuller[20b8| ), sub-fragmentation 

of prestellar cores into binary or multiple systems complicates the di¬ 
rect mapping of the prestellar CMF onto the IMF of individual stars. 
Lacking sufficient spatial resolution to probe core multiplicity with the 
present Herschel observations, we do not enter this debate here and con¬ 
centrate on the relationship between the prestellar CMF and the system 
IMF. 
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Mass, M (Me) 


Fig. 17. Comparison of the CMF observed for the majority (~81%) 
of candidate prestellar cores lying within the 0.2 pc-wide footprints of 
the DisPerSE-detected filaments (light blue histogram) to that observed 
for the minority ~19%) of prestellar cores lying outside these filaments 
(magenta histogram). The upper dark blue histogram and the other lines 
are the same as in Fig.p^ 


prestellar CMF observed off filaments may peak at a somewhat 
lower mass. A two-sample K-S test indicates that the probability 
that the CMFs observed on and off filaments above 0.2 M© are 
drawn from the same intrinsic distribution is only ~ 2% (equiv¬ 
alent to a ~ 2.3 (T result in Gaussian statistics). This is not a 
very strong conclusion, however. First, there are only 54 can¬ 
didate prestellar cores with masses > 0.2 M© lying outside the 
0.2 pc-wide filament footprints, implying that our estimate of the 
prestellar CMF off filaments suffers from small-number statis¬ 
tics. In fact, we cannot even exclude the possibility that some 
of the prestellar cores presently classified as lying off filaments 
may be associated with faint filaments not identified with Dis- 
PerSE in Sect. |4.2[ Second, the median background cloud col¬ 
umn density observed off filaments is lower (A^^^ ~ 4) than the 
median background cloud column density observed on filaments 
(^back _ y Accordingly, the completeness level of our Her- 
schel survey for prestellar cores is expected to be somewhat bet¬ 
ter off filaments than on filaments (see Fig. |B.6| in Appendix B), 
which may slightly bias the direct comparison of the two CMFs. 


5.6. A quasi-universal efficiency of the star formation 
process in dense gas? 

Our Herschel results on the prestellar core formation efficiency 
(CFE) as a function of column density in the Aquila cloud (see 
Sect. 5.2 and Fig.p^ connect very well with recent near-/mid- 
infrared studies of the star formation rate (SFR) as a function 


of gas surface density in nearby molecular clouds (e.g., Heider- 
|man et aL]|2010t |Lada et J?||2010[ |2012t [Evans et al.||2014| ). 
These infrared studies show that the global SFR derived from 
direct YSO counting (as opposed to the prestellar core count¬ 
ing used in the present study) tends to be linearly proportional 
to the mass of dense gas above a surface density threshold 
corresponding to ~7-8, and drops to much lower val¬ 

ues below the threshold. This column density threshold is es¬ 
sentially the same as that found with Herschel for the forma¬ 
tion of prestellar cores in the Aquila cloud (cf. Figs. m &[g. 
Moreover, the star formation rate per unit mass of dense gas 
above the threshold found by infrared studies of nearby clouds. 

Article number, page 22 of[^ 


yr dLada et al.l2010l|2012| 


namely SFR/Mdense ~ 4.6 x 10“ 

or SFR/Mdense ~ 2.5^|-^ X 10“^yr“^ ([Evans et al.||2014|), is en 
tirely consistent with the roughly constant prestellar CJFE de¬ 
rived for Av > 7 in Aquila, which corresponds to SFR/Mdense ~ 
5^2 X 10“^ yr“^ (see horizontal dotted line in Fig.f^ adopting a 
typical prestellar core lifetime = 1.2 Myr (see 5ect. |5.1| ) and 
a loc al st ar formation efficiency 6core = 0-4 at the core level (see 
Sect.[^. _ __ 

As pointed out by [Lada et al.| ( |201()||2012| ), the nearby cloud 
value of the “efficiency” of the star formation process in dense 
gas is also very similar to the efficiency value S FR/Mdense ~ 
2 X 10“^yr“^ found by Gao & Solomon ( |2Q04 ) for external 
galaxies, using HCN observations of dense gas and far-infrared 
(IRAS) estimates of the SFR in galaxies. While direct compar¬ 
ison between the Galactic and extragalactic values is affected 
by large uncertainties because different tracers of dense gas and 
star formation were used by [Lada et al.| ( |201Q| ) on the one hand 
and [Gao & Solomonj ( |2004| ) on the other, these results suggest 
that there may be a quasi-universal “star formation law” within 
dense gas above the (column) density threshold. Equivalently, in 
terms of a concept often used in the extragalactic community, 
this means that there may be a quasi-universal depletion time, 
= Mdense/SFR ~ 20-50 Myr, for the dense gas above the 
threshold. This “star formation law” is not strictly universal since 
its does not seem to apply to the extreme environmental condi¬ 
tions of the central molecular zone near the Galactic center, for 
instance, where star formation is observed to be more inefficient 
above the density threshold, by more than an order of magnitude 
( jLongmore et al.|2()T3] ). 

Our Herschel findings in the Aquila cloud allow us to go 
one step further and link this quasi-universal efficiency of the 
star formation process in dense gas to three parameters charac¬ 
terizing the physics of prestellar cores, i.e., the core formation 
efficiency in supercritical filaments, /pre, the lifetime of prestel¬ 
lar cores, ^pre, and the efficiency of the conversion from prestellar 
core mass to stellar system mass, 6core, i-e., 


(see jAndre et al.|20T^ . 


6. Summary and conclusions 

We used the SPIRE and PACS parallel-mode maps taken as part 
of the Herschel Gould Belt survey to obtain an extensive cen¬ 
sus of dense cores and their connection with molecular cloud 
structure in the Aquila star-forming region. Our main results and 
conclusions may be summarized as follows: 

1. The high-resolution (~ 18" or ~ 0.02 pc) column density 
map that we derived from the Herschel photometric data 
shows that the Aquila cloud is highly filamentary and fea¬ 
tures a column density probability density function (PDF) 
with a prominent power-law tail above Ay ~ 5-7. About 
10 %-20% of the gas mass is in the form of filaments below 
Ay ~ 7, while as much as ~ 50%-75% of the gas mass is in 
the form of filamentary structures above Ay ~ 7-10. 

2. In the ~11 deg^ field imaged with both SPIRE and PACS at 
five wavelengths from 70 pm to 500 yum, we identified 651 
starless cores, 446 candidate and 292 robust prestellar cores, 
and 58 protostellar cores (such as Class 0 objects), based 
on multi-scale, multi-wavelength core extraction with the 
getsources algorithm. The samples of candidate and robust 
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prestellar cores were estimated to be ~ 90% and ~ 80% com¬ 
plete, respectively, down to an observed core mass ~ 0.2 M©. 
The candidate prestellar cores have estimated median mass 
~ 0.45 Mq, median deconvolved FWHM diameter ~ 0.03 pc, 
median average column density ~ 7 x 10^^ cm“^ and median 
average volume density ~ 4 x lO'^ cm“^. 

3. The typical lifetime of the Herschel prestellar cores was esti¬ 
mated to be ^pre = 1.2+ 0.3 Myr or ~ 4 free-fall times (^ff) and 
to decrease from ^pre -1.4 Myr for cores with average vol¬ 
ume density ^10^ cm“^ to a few times 10"^ yr for cores with 
average volume density ^10^ cm“^. The densest prestellar 
cores in the sample appear to have a lifetime comparable to 
their free-fall timescale and may be collapsing. 

4. There is strong evidence of a column density threshold for 
the formation of prestellar cores, at an equivalent visual ex¬ 
tinction level ~ 7, in the sense that the probability func¬ 
tion of finding a prestellar core increases by more than an 
order of magnitude from ~ 4 to ~ 10 and is well fit 
by a smooth exponential step function. Likewise, the prestel¬ 
lar core formation efficiency (CFE) or fraction of cloud mass 
in the form of prestellar cores was found to increase by about 
two orders of magnitude between A^ ~ 5 and A^ ~ 15 and 
to reach a roughly constant value CFEmax = /pre ~ 15% at 
higher column densities. This reflects a significantly sharper 
transition than predicted by “multi-freefall” models of the 
star formation rate in molecular clouds, and argues for the 
presence of a true physical (column) density threshold for 
prestellar core formation. 

5. The compact dense cores are closely associated with the fila¬ 

mentary structure, and preferentially the densest filaments. In 
particular, a very high fraction (75%^5^^^'’) of prestellar cores 
were found to lie within supercritical filaments with masses 
per unit length Mune > Miine,crit, where Miine,crit = 2c^/G ~ 
16 Mo/pc is the critical mass per unit le ngth of nearly isother- 
mal, long cylinders at T ~ 10 K (see llnutsuka & Miyamal 
[T997] ). - 

6 . The prestellar core mass function (CMF) derived using the 
samples of 446 candidate and 292 robust prestellar cores is 
well fit by a lognormal distribution, peaks at ~ 0.4-0.6 M©, 
and is very similar in shape to the system IMF. This CMF is 
consistent with an essentially one-to-one mapping between 
prestellar core mass and stellar system mass with a local 
star formation efficiency 6core ~ 0-4^o i within an individual 
prestellar core. 

7. Our Herschel findings in the Aquila cloud connect very 
well with recent Spitzer studies of the star formation rate in 
nearby molecular clouds. They support the view that there 
may be a quasi-universal “efficiency” of the star formation 
process in dense gas, SFR/Mdense ~ X 10“^yr“^ and 
that this quasi-universal “efficiency” may be closely linked 
to the physics of prestellar core formation within filaments: 

SFI^/Afdense ~ ./pre X ^core/^pre* 
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Appendix A: A catalog of dense cores identified 
with Herschel in the Aquila cloud complex 

Based on our Herschel SPIRE/PACS parallel-mode imaging sur¬ 
vey of the Aquila cloud complex, we identified a total of 749 
dense cores, including 685 starless cores and 64 protostellar 
cores. (Among these, 34 starless cores shown as yellow triangles 
in Fig.[^ as well as 6 protostellar cores, were excluded from the 
scientific discussion of Sect. due to likely contamination by 
more distant, background objects - see Sect. |4.3[ ) The master 
catalog listing the observed properties of all of these Herschel 
cores is available in online Table |A.l] A template of this online 
catalog is provided below to illustrate its form and content. 

The derived properties (physical radius, mass, SED dust tem¬ 
perature, peak column density at the resolution of the 500 //m 
data, average column density, pe ak vo lume density, and average 
density) are given in online Table |A.2] for each core. A portion of 
this online table is also provided below. The derived properties 
of the Herschel-doitci^d protostars and YSOs will be published 
in a forthcoming paper. 
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Table A.l. Catalog of dense cores identified in the RGBS maps of the Aquila complex (template, full catalog only provided online). 
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run No 303 


run No 79 


Aquila core HGBS^ 183005.4-014833 core HGBS^ 182754.3-034237 

Tdurt (K) - 9.2±0.6. Moss (Mg) - 1.01 ±0.28 (K) = 12.3±0.5, Mass (Mg) = 2.0610.26 




Wavelength (/xm) Wavelength {firry) 


Fig. A.l. Examples of Herschel spectr al en ergy distributions (SEDs) for a prestellar core (left, see Eig. |A.3| for the corresponding image cutouts) 
and a protostellar core (right, see Eig. |A.4| for the corresponding image cutouts). These SEDs are constructed from the background-subtracted 
integrated flux densities (cross symbols) measured by getsources. A graybody fit to the SED observed longward of 160 pm is superimposed as a 
blue curve in both panels. Only upper limits are available at 70 pm and 160 pm for the prestellar core shown in the left panel. Similar SED plots 
are provided online for all selected cores. A single-temperature graybody rarely provides a good fit to the overall SED of a protostellar core but 
can nevertheless describe the SED longward of 160 pm reasonably well (cf. right panel). 
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Fig. A.2. Blow-up column density images of two Aquila subfields at 18.2" resolution. Black and red ellipses mark the EWHM sizes of the starless 
cores and protostellar cores, respectively, sel ected from getsources e xtractions in these two subfields. Green ellipses show the EWHM sizes of the 
sources independently detected with CSAR ( |Kirk, J. M. et al.|2013| t in the high-resolution column density image. 
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Aguila core HGBS^ 103005.4-0146^3 
30 05,46 Dm —01 49 33,6 



Aquila core HGBS_J 185005.■*1-014833 
Rfl = l0 50 05.46 Dec=-01 43 55.& 



run Nq 303 

Aqijilt] core HGBSJ 18.3005.4-0 1493.5 
Ra-t8 30 05.46 Dec—01 48 35,6 



Aquila core HGeS_Ji85005.4-014335 
Ra=lS 50 05.46 Dec=-0l 40 33.6 




Fig. A.3. Example blow-up Herschel images at 70/160/250/350/500 jjim and high-resolution column density map for a (bound) prestellar core. 
Ellipses represent the estimated major and minor EWHM sizes of the core at each wavelength; they are show n as solid or dashed curves depending 
on whether the core is significantly detected or not, respectively, at a particular wavelength. See Table |A.2] for the physical radius of the core and 
other derived properties. Similar image cutouts are provided online for all selected starless cores. 


run No 79 



Aquila cors HGBS_JlS2754.3-034237 

Ra=ie 21 54.36 Dec=-03 42 37.2 


Aquila corfl HCB£_J 182754.3-034257 
Ra=l8 27 54.56 Dec=-03 42 57.2 
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Fig. A.4. Same as Eig.|A.3|for a protostellar core. Similar image cutouts are provided online for all selected protostellar cores. 
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Appendix B: Completeness of HGBS prestellar core 
extractions in Aquila 

To estimate the completeness of our census of prestellar cores in 
Aquila, we used several sets of simulated data on the one hand 
(Sect. 113, and a simple model of the core extraction process 
and completeness problem on the other (Sect. El- 



Fig. B.l. Synthetic core mass function (CMF) derived from simulated 
source extractions (blue histogram) compared to the input mass function 
(dA/dlogM oc of a population of 5622 model cores (red curve) 

constructed as described in the text. The estimated 90% completeness 
level (in observed core mass) is indicated by the vertical dashed line at 
0.2 Mq. The black solid line shows a power-law fit to the derived CMF 
above the 0.2 M© completeness level; it is in excellent agreement with 
the input core mass function. The drop of the synthetic CMF below the 
input CMF at the high-mass end is due to the fact that the derived core 
masses t end t o underestimate the true core masses by ~ 20% on average 
(see Fig. 


Appendix B.1: Monte-Carlo simulations 


To simulate real core extractions, we first constructed clean maps 
of the background emission at all Herschel wavelengths (includ¬ 
ing a column density plane), by subtracting the emission of the 
compact c ores id entifi ed with getsources in the observed data 
(cf. Sects. |4.4| & |4.5| ). We then inserted several sets of model 
Bonnor-Ebert-like cores throughout the clean-background im¬ 
ages in order to generate a full set of synthetic Herschel and 
column density images of the region. In the example illus¬ 
trated in Figs. |B.1| & |B.2[ for instance, we used a population 
of 5622 model starless cores with a flat input mass distribu¬ 
tion (dA/dlogM oc similar to the mass distribution of CO 

clumps) from 0.02 M© to ~ 30 Mq. This example is particularly 
useful as it allowed us to test the robustness of the conclusion 
that the observed prestellar CMF is signiflcantly steeper than the 
mass distribution of CO clumps. The model cores had positions 
in a mass versus size diagram consistent with critical Bonnor- 
Ebert isothermal spheres at effective gas temperatures ~ 7-20 K. 
The dust continuum emission from the synthetic Bonnor-Ebert 
cores in all Herschel bands was simulated using an extensive 
grid of spherical dust radiative transfer models constructed by 
us with the MODUST code (e.g., Bouwman et al.|20()m|Bouw-] 
|man||200f] ). In particular, each of the synthetic prestellar cores 
was given a realistic dust temperature profile with a significant 
drop in dust temperature toward core center, as observed in the 
case of spatially-resolved starless cores (cf. jRoy et al.|2014| ). The 
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synthetic cores were spatially distributed randomly over the re¬ 
gions of the column density map where > 5 x 10^^ cm“^ 
(containing most, if not all, of the observed prestellar cores in 
the real data - see Sect. |53, with no particular mass segrega¬ 
tion. Once satisfactory synthetic skies resembling the observed 
images had been generated, compact source extraction and core 
selection/classiflcation were performed with getsources in the 
same way as for the real data (see Sects. [4^&[45l). 

As mentioned in Sect. |4.8| and shown in Fig. the results 
of these Monte-Carlo simulations suggest that our Herschel cen¬ 
sus of prestellar cores in the Aquila cloud complex is ~ 90% 
complete down to ~ 0.3 M© in true core mass. Figure |BT] fur¬ 
ther illustrates that the core mass function can be reliably de¬ 
termined down to the completeness mass limit. In this example, 
a Kolmogorov-Smirnov (K-S) test shows that the derived CMF 
is statistically indistinguishable (at the ~ 90% confidence level) 
from the input mass function above the completeness limit. In 
particular, the best-fit power -law function to the derived CMF 
(black solid line in Fig. 113 is identical to the input dA/dlogM 
oc M~^-^ power law. This test therefore confirms that the best- 
fit power law to the observed CMF (dA/dlogM oc -33+0.06 _ 
see Sect. |5.5| ) is signiflcantly steeper than the typical mass dis¬ 
tribution of CO clumps/clouds (dA/dlogM oc - e.g. Blitz 

1993; Kramer et al. 1998) and cannot be an artifact of the core 
extraction process. 

The same Monte-Carlo simulations were also used to assess 
the accuracy of the main derived parameters (e.g. core mass, ra¬ 
dius, and dust temperature) by comparing the estimated values 
after core ext racti on to the intrinsic input values of the model 
cores. Figure |B.2^ shows that the derived core masses tend to 
underestimate the true core masses by -20-30% on average, 
and Fig. shows that the derived SED temperatures tend 
to overestimate the intrinsic mass-averaged dust temperatures of 
the cores by typically ~ 1 K. A similar plot for the core sizes 
(Fig.|B3 suggests that the derived core sizes (prior to decon¬ 
volution) are quite reliable and remain within - 5% of the true 
core sizes on average. We interpret the mass effect ( Fig.|B.2^ ) as 
a direct consequence of the temperature effect (Fig. |B.2} )) since 
overestimating the dust temperatures leads to underestimating 
the core masses. The temperature effect arises from the fact that 
the dust temperature derived from a global fit to the SED of a 
starless core overestimates the mass-averaged dust temperature 
owing to a distribution of dust temperatures along the line of 
sight (see |Roy et al.|2014| and Sect. |4.6| ). 

Taking the ~ 20-30% mass effect into account, we conclude 
that the - 90% completeness limit at ~ 0.3 Mq in true core mass 
corresponds to - 0.2 Mq in observed core mass. 


Appendix B.2: Model of the completeness problem 


The Monte-Carlo simulations described above provide an esti¬ 
mate of the global completeness limit of the core survey. The 
completeness level of the core extractions is, however, expected 
to be background dependent. To assess the importance of this 
dependence, we constructed a simplified model of the core ex¬ 
traction process. 

Owing to the high sensitivity and quality of the Herschel im¬ 
ages, the HGBS survey is not limited by instrumental noise but 
by confusion arising from small-scale cloud structure, an effect 
commonly referred to as “cirrus confusion noise” in the literature 
(e.g., see |Gautier et al.|19^|Kiss et al.|2Q01t|Roy et al.|2010 ). 
To estimate the level of such cirrus confusion noise from the 
Herschel data, we measured the rms level of background fluctu- 
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Fig. B.2. Left: Ratio of derived to intrinsic (or ‘true’) core mass as a function of derived core mass for the same set of simulated core extractions 
as used in Sect. |4.8| and Fig. |B.l| The error bars are +lcr where cr is the dispersion of the mass ratio in each mass bin. The median mass ratio is 
~ 0.8 above 0.4 M© (as indicated by the horizontal blue line) and ~ 0.7 close to the 90% completeness limit of 0.2 M© in observed core mass. The 
horizontal dashed line marks the mass ratio of 1 expected in the case of perfect core extractions and mass estimates. Right: Difference between 
derived SED temperature and intrinsic mass-averaged dust temperature as a function of derived core mass for the same set of simulated core 
extractions. The error bars are +lcr where cr is the dispersion of the temperature difference in each mass bin. The median temperature difference is 
about -hO.8 K above 0.4 M© (as indicated by the horizontal blue line) and ~ 1 K close to a derived core mass of 0.2 M© (completeness limit). The 
horizontal dashed line marks the zero difference expected in the case of perfect core extractions and temperature estimates. 



Fig. B.3. Ratio of derived to true core size as a function of deri ved c ore 
mass for the same set of simulated core extractions as in Fig. |B.l| and 
Fig. |B.2| The error bars are +lcr where cr is the dispersion of the size ra¬ 
tio in each mass bin. The horizontal dashed line marks the size ratio of 1 
expected in the case of perfect core extractions and size estimates. Note 
how the median core size measured in each mass bin remains within 
5% of the true core size above the ~ 90% completeness of ~ 0.2 M© in 
derived core mass. 


ations in a sliding box V x V in siz over the entire column 
density map of the Aquila complex after subtracting the sources 
identified by getsources. Correlating the resulting map of rms 
fluctuations with the input background column density map led 
to Fig. B.4 which clearly shows that the level of column density 


The size of the sliding box corresponds to ~ 0.075 pc x 0.075 pc at 
d ~ 260 pc, which is similar to the size scale of prestellar cores. 


fluctuations increases with background column density approxi¬ 
mately as a power law: 


A^H2,nns~ 3.9 X 10^0 cm-^x 


^Hz.back 

7 X 10^^ cm“^ 


1.6 


(B.l) 


The power-law index of 1.6 derived here from Herschel data is 
very similar to that reported in earlier papers discussing cirrus 
noise (e.g. Gautier et al. 1992, Kiss et al. 2001, Roy et al. 2010). 
Since the level of background fluctuations increases with column 
density, one expects core extraction to be increasingly more dif¬ 
ficult and thus survey completeness to decrease significantly in 
higher column density areas within the field. 

The model we used to estimate the magnitude of this effect 
and get around the problem of a background-dependent com¬ 
pleteness level was based on the following assumptions: 

• A dense core is defined as the immediate vicinity of a 
column density peak departing significantly, i.e., by more than 
5 X A^H 2 ,rms the field of background cloud fluctuations (see 
core selection criteria in Sect |4.5| ). 

• A prestellar core, i.e., a self-gravitating starless core, can 
be approximately modeled as a critical Bonnor-Ebert spheroid 
of mass Mbe and outer radius Rbe, bounded by the grav itational 
pressure of the background cloud Rback 


0.88 G dMcKee 


& Tan|[2003] ), where Sback = x A^H 2 ,back- The mean mtrm 


sic column density contrast of such a model prestellar core is 
^BE/^back ~ 1-5, where Ebe = ^be/(^^be^' 

• The ability to detect a core in the Herschel data depends 
primarily on the apparent column density significance of the 
core defined as Score,obs/^rms, where Score,obs is the apparent (ob¬ 
served) column density of the core after convolution with the 


observing beam, i.e., S, 


'core,obs 


— ^core / 


core,conv 


), and Znns = 


are 


//mHXA^H 2 ,rms- The Monte-Carlo simulations of Appendix B.l 
consistent with this assumption and suggest that the complete 
ness level is ^ 90% for cores with an apparent column density 
significance larger than 5 (see Fig. |B3J. 
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Background cloud column density (H^cm 



Apparent significance over background fluctuations 


Fig. B.4. Median root mean square (rms) level of background column 
density fluctuations as a function of background cloud column den¬ 
sity as measured in a 1' x 1' sliding box over the clean background 
image of the Aquila complex produced by getsources from the high- 
resolution column density map. The error bars correspond to the in¬ 
terquartile range of background fluctuations about the median level in 
each column density bin. The straight line represents a weighted power- 
law fit to the data points above N 112 ~ 2 x 10^^ cm“^. 


In outline, our simplified model of the completeness problem 
may be described as follows: 

• Two effects, beam dilution and temperature dilution, can 
make the apparent column density contrast SBE,obs/^back of a 
model core smaller than its intrinsic column density contrast of 
1.5: 

^BE,obs/^back ~ 1.5 X {RbE / ^BE,com) ^ [^Veff(^core)/^yeff(^back)? 

where Veff is a fiducial Herschel observing frequency which 
we take to correspond to T ~ 350yum. Taking advantage of 
the fact that the column density distribution of a Bonnor-Ebert 
core with outer radius Rbe is well approximated by a Gaus¬ 
sian distribution of FWHM ~ Rbe, the observed radius of the 

core is approximately RBE,com = (R^e HPBW (where 
HPBW corresponds to the half-power beam width resolution of 
the column density map projected at the distance of the Aquila 
cloud), and the beam dilution factor can thus be expressed as 

(/?BE//?BE,conv)2 = 1/[1 + (HPM/RseH 

• The apparent column density significance can be written as 
the product of the apparent column density contrast and a cirrus 
noise factor, SBE,obs/^rms ~ (^BE,obs/^back) ^ (^back/^rms)? whcrC 
the cirrus noise factor is: 


^back/^rms ~ ^^H2,back/A^H2 


18 X 


A^H2,back 
7 X 10^^ cm“^ 


-0.6 


according to Eq. (B.l). 

• Assuming that the fundamental completeness curve is the 
completeness function 7^( 5) o f apparent column density sig¬ 


nificance S shown in Eig. B.5 completeness can be estimated 


as a function of core mass and background column density as 
C(MBE,Sback) = ^r[5(MBE,Sback)]. The Corresponding function 
of Mbe is shown for five values of the background column den¬ 
sity A^H 2 ,back in Eig. |B.6[ Figure [B^ shows how the completeness 
of prestellar core extractions is expected to decrease as back¬ 
ground cloud column density and cirrus noise increase. 
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Fig. B.5. Completeness curve as a function of apparent column density 
significance over local background cloud fluctuations derived from the 
Monte-Carlo simulations described in Sect. 14.81 

Model completeness curves at d = 260pc (Aquila) 



True core moss (Mq) 


Fig. B.6. Model completeness curves of Herschel prestellar core extrac¬ 
tions in Aquila for five values of the background cloud column density 
expressed in units of visual extinction from Ay,back = 5 to Ay back = 20. 


• To estimate a global completeness curve for our census 
of prestellar cores in the Aquila complex, we used the observed 
distribution of mass in the cloud as a function of background 
column density (cf. Figs, [^/b) and took advantage of the ex¬ 
istence of a column density “threshold” at Ay,back ~ 5-7, above 
which the bulk of core and star formation is believed to occur (cf. 
Sect. |5.2| and Fig. ED and the column density PDF is well fitted 
by a power-law distribution. We also assumed that the number 
of prestellar cores in the cloud scales linearly with cloud mass 
above the threshold. This assumption is consistent with recent 
infrared studies which find that the global star formation rate 
tends to be linearly proportional to the mass of dense gas above 
the threshold (e.g.,[Heiderman et al.|20rm|Fada et al.|2010t|Gao 

& Solomon||2004| ). It is also consistent with the roughly con¬ 

stant prestellar core formation efficiency found here above the 
threshold (see Fig.[^. The global completeness curve was thus 
computed as a weighted average of the individual completeness 
curves at fixed background column densities: 


@C{Mbe) = 


^^dei 


/ 


C(Mbe, S^back)-(Av,back)^Av,baek- 

'Av=5 
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The resulting global completeness curve, which represents 
the best estimate of the completeness of our Herschel survey 
for prestellar cores in Aquila according to our model, is shown 
in Fig. It can be seen that this global completeness curve 
is very similar to the individual completeness curves for back¬ 
ground column densities close to the threshold (see Ay, back = 5- 
10 curves in Fig. |B.6| ). It is also very similar to the empirical 
completeness curve derived from Monte-Carlo simulations (see 
Sect. |4.8| ). The model completeness curve is almost flat above a 
true core mass level of 0.3 M©. Using this model curve to cor¬ 
rect the observed CMFs of candidate and robust prestellar cores 


for incompleteness would only have a minimal effect in Fig. 16 
above an observed core mass level of ~ 0.2 M©. (The corrected 
CMFs differ from the uncorrected CMFs only below ~ 1 M© and 
by much less than the uncertainty area displayed in light blue in 

Fig. fra) 


Global completeness curve at d = 260 pc (Aquila) 



True core mass (M©) 


391, while the total number of starless cores (651) would re¬ 
main the same. Accordingly, the estimated lifetime of candidate 
prestellar cores would also slightly increase from ~ 1.4 Myr to 
~ 1.8 Myr, and that of robus t prestellar cores from ~ 0.9 Myr 
to ~ 1.3 Myr (see Sect. 15.11), leading to t^^Q = 1.5 + 0.3 Myr. 
The prestellar core formation efficiency (CFE) as a function of 
background column density (cf. Fig. [T^, and in particular the 
roughly constant value CFEmax = /pre ~ 15% at high col¬ 
umn densities, would not change. Our corresponding estimate 
of th e “ef ficiency” of the star formation process in dense gas (cf. 
Sect. 5.6), SFR/Mdense = /preX^core/^pre, would tiowever decrease 
from 5 X 10“^ yr“^ to 2 x 10“^ yr“\ becoming closer to the effi¬ 
ciency value reported by |Evans et al.|(|2Q14|) and|Gao & Solomon 


Finally, the 
T3l and 


( |2004| ) than to the value found by |Lada et al.] ( [26l0|) . Fin 
column density maps shown in Fig.[lJ Fig.|^ Fig.|^ Fig. 

Fig.[^ as well as the spatial correspondence between cores and 
filaments, would remain unchanged. The scaling of our column 
density maps in terms of mass per unit length along the filaments 
would however change by ~ 60% upward, since the character¬ 
istic physical width of the filaments would increase by ~ 60%. 
As a consequence, the white areas which highlight supercritical 
filaments in Figs. and would slightly expand, im¬ 

proving the correspondence between the spatial distribution of 
prestellar cores/protostars and that of supercritical filaments. To 
summarize, our main conclusions do not depend strongly on the 
adopted distance. 


Fig. B.7. Model completeness curve of Herschel prestellar core extrac¬ 
tions in the Aquila cloud complex {d - 260 pc) as a function of intrinsic 
model core mass. 


Appendix C: Effect of distance uncertainty 

As mentioned in Sect. there is some ambiguity concerning 
the distance to the Aquila molecular cloud comple x. A number 
of arguments, presented by |Bontemps et al.| ( |2010| ) and summa¬ 
rized in Sect. suggest that the bulk of the region studied here 
and shown in Fig. corr esponds to a coherent c loud complex 
at d- = 260 pc (see also |Gutermuth et al.|[2008| ), which is the 
default distance adopted in the present paper. Other studies in 
the literature (see references in Sect.j^, however, place the com¬ 
plex at the larger distance, d+ = 415 pc, of the Serpens Main 
cloud ( |Dzib et al.|2010 ). It is thus worth discussing how our re¬ 
sults would be affected if we had adopted the larger distance esti¬ 
mate, J+, instead of J_. The core mass estimates, which scale as 
Syd?-l{By{T^)Ky\ where Sy is integrated flux density and By{T^) 
is the Planck function, would systematically increase by a factor 
of 2.5. This would shift the CMFs shown in Fig.fraand Fig.pT 


to the right and thus lower the efficiency 6core from 0.4^^^ to 
0.2 + 0.1. In comparison, the core size estimates, which scale 
linearly with distance d, would increase by only 60%. The BE 
mass ratio ckbe = AfBE,crit/Afobs 5 listed in Col. 17 of online Ta¬ 
ble A.2, scales as d~^ and would decrease by 60% for all cores. 
Accordingly, all cores would move upward as indicated by an 
arrow in the mass versus size diagram of Eig.[^ which would in¬ 
crease the fraction of prestellar cores among starless cores from 
60% + 10% to 70% + 10%. More precisely, the number of can¬ 
didate prestellar cores would increase from 446 to 565 and the 
number of robust prestellar cores would increase from 292 to 
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Fig. 18. Distribution of mean FWHM inner widths for the 90 filaments 
traced with DisPerSE in the Aquila entire field (see blue skeleton in 
Figs.|^&[^and Sect. |4.2| ). These widths results from a filame nt profile 
analysis similar to that described in Arzoumanian et aL| ( |201 1| ) and were 
deconvolved from the 18.2" HPBW resolution of the high-resolution 
column density map used to construct the radial profiles of the filaments. 
The median filament width is 0.12 pc, as marked by the vertical dotted 
line, and the standard deviation of the distribution is 0.04 pc. 
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Fig. 19. a) Distribution of dust temperature values in the Aquila temperature map shown in Fig. |b) Distribution of SED dust temperatures for 
all selected starless cores with reliable SED fits (see Sect. |4.6| ). Note how the distribution of core temperatures peaks at significantly lower values 
than the distribution of background cloud temperatures. 





















